Mon. Not. R. Astron. Soc. 000, ITH291 f????1 



Printed 31 May 2011 



(MN IATeX style file v2.2) 



The Proto- Magnet ar Model for Gamma-Ray Bursts 



B. D. Metzger 1 ' 5 *, D. Giannios 1 , T. A. Thompson 2 , N. Bucciantini 3 , E. Quataert 4 

1 Department of Astrophysical Sciences, Peyton Hall, Princeton University, Princeton, NJ 08544> USA 

2 Department of Astronomy and Center for Cosmology & Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA 
3 NORDITA, AlbaNova University Center, Roslagstullsbacken 23, SE 10691 Stockholm, Sweden 

4 Astronomy Department and Theoretical Astrophysics Center, 601 Campbell Hall, University of California, Berkeley, CA 94720, USA 
5 NASA Einstein Fellow 



Accepted . Received ; in original form 



ABSTRACT 

Long duration Gamma-Ray Bursts (GRBs) originate from the core collapse of mas- 
sive stars, but the identity of the central engine remains elusive. Previous work has 
shown that rapidly spinning, strongly magnetized proto-neutron stars ('millisecond 
proto-magnetars') produce outflows with energies, timescales, and magnetizations cro 
(maximum Lorentz factor) that are consistent with those required to produce long du- 
ration GRBs. Here we extend this work in order to construct a self-consistent model 
that directly connects the properties of the central engine to the observed prompt 
emission. Just after the launch of the supernova shock, a wind heated by neutrinos 
is driven from the proto-magnetar. The outflow is collimated into a bipolar jet by 
its interaction with the progenitor star. As the magnetar cools, the wind becomes 
ultra- relativistic and Poynting-flux dominated (oo 3> 1) on a timescale comparable 
to that required for the jet to clear a cavity through the star. Although the site and 
mechanism of the prompt emission are debated, we calculate the emission predicted 
by two models: magnetic dissipation and shocks. 

Magnetic reconnection may occur near the photosphere if the outflow develops 
an alternating field structure due to e.g. magnetic instabilities or a misalignment be- 
tween the magnetic and rotation axes. Shocks may occur at larger radii because the 
Lorentz factor of the wind increases with time, such that the faster jet at late times 
collides with slower material released earlier. Our results favor magnetic dissipation 
as the prompt emission mechanism, in part because it predicts a relatively constant 
'Band' spectral peak energy -E pca k with time during the GRB. The baryon loading 
of the jet decreases abruptly when the neutron star becomes transparent to neutri- 
nos at t = t„-thin ~ 10 — 100 seconds. Jets with ultra-high magnetization cannot 
effectively accelerate and dissipate their energy, which suggests this transition ends 
the prompt emission. This correspondence may explain both the typical durations of 
long GRBs and the steep decay phase that follows. Residual rotational or magnetic 
energy may continue to power late time flaring or afterglow emission, such as the 
X-ray plateau. We quantify the emission predicted from proto-magnetars with a wide 
range of physical properties (initial rotation period, surface dipole field strength, and 
magnetic obliquity) and assess a variety of phenomena potentially related to magnetar 
birth, including low luminosity GRBs, very luminous GRBs, thermal-rich GRBs/X-ray 
Flashes, very luminous supernovae, and short duration GRBs with extended emission. 

Key words: Stars: neutron; stars: winds, outflows; gamma rays: bursts; MHD 



1 INTRODUCTION 

Soon following the disc overy of Gamma-Ray Bursts (GRBs; 
iKlebesadel et all 1 19731 ). there were possibly mo re theories 
for their origin than theorists (jRudermanl Il97a ) . However, 

* E-mail: bmetzger@astro.princeton.edu 



once GRBs w ere confirmed to orig inate from cosmological 
distances (e.g. iMetzger et all 1 19971 ) ■ the joint requirements 
of supernova-scale energies, short (millisecond) timescales, 
and relativistic speeds significantly narrowed the list of plau- 
sible central engines. It is now generally accepted that GRBs 
result from the formation or catastrophic rearrangement of 
stellar- mass black holes (BHs) or neutron stars (NSs). This 
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Figure 1. Schematic diagram of the regimes of neutron star versus black hole formation in core collapse SNe at sub-solar metallicities 
(solid line) in the space of main sequence mass and initial proto-NS spin period Ppj, taking into account the possible effects of rapid 
rotation and strong magnetic fields. The dotted line denotes the rotation rate above which the NS rotational energy E TO t (eq. [I]) exceeds 
the gravitational binding energy of the progenitor envelope. The dashed line denotes the rotational energy E I0 t = 10 52 ergs sufficient to 
power a 'hypernova'. The right axis shows the magnetic field strength B^ip that would be generated if the magnetic energy in the dipole 
field is ~ 0.1% of E TO t (eq. [?])• The dot-dashed line is the minimum rotation rate required for a magnetar with a field strength i?di P to 
produce a classical GRB with energy E 1 > 10 51 ergs, based on the model presented in [|4] 



conclusion has only been strengthened in recent years due to 
the much richer picture of the prompt and afterglow emis- 
sion provided by the Swift and Fermi missions. However, 
despite a wealth of new data, the identity of the central 
engine remains elusive. 

At least some long duration GRBs originate from th e 
deaths of very massive stars ()Wooslev fc Blooml [2006), 
as confirmed by their observed associ ation with energetic 



Galama et all 1998: 



Chornock et al.l 



core collapse super n ovae (SNe) I 

Bloom et all 1 19991: IStanek et al.l 12003 
20ld : IStarling et all l2010r ). It nevertheless remains 

unsettled wheth er the c entra l engine is a rapidly 

accre ti ng BH JWooslevi Il993l: iMacFadven fc Woosle 



1999 ; iNagataki et all 120071 ; iBarkov fc Komissarovl 12*00 
Lindner et all l2010l ) or a rapidly spinni ng, strongly 
magnetized NS ( a 'millisecond magneta r' ; lUsovl Il992 | ; 
Thompsonl Il994l; iBlackman fc Yjl Il998l; IWheeler et all 



200d: I Zhan g feMeszarosI l200ll; iThompson et all I200J 



Metzger et 



ng fe Me 
all 120071 : 



Bucciantini et all 120071 . 120081 . l2009h . 



Although much less is known about the origin of short 



duration GRBs, the properties of their host galaxies and 
their notable lack of an accompanying SN are consistent 
with an origin ass ociated with the merger of NS-NS and 
NS-BH binaries dHiorth et all 120051: iBloom et all 120061 : 
iBerger et al.ll2005l ; see e.g. Bergerll2010l for a recent review). 
However, the unexpected discovery that many short GRBs 
are followed by an energetic X-ray 'tail' lasting ~ 100 sec- 
onds has challenged bas i c predictions of the merger model 
(e.g. iGehrels et al.l [20061 : iGal-Yam et~all 120061 : IPerlev et all 
l2009h and may hint at an alternative origin for some events, 
such as magnetar formation via the accretion-in duced 
collapse (AIC) of a white dwarf ijMetzger et alll2008l ). 

The large range in length scales and the complexity of 
the physics involved in producing a GRB have thus far pre- 
vented all steps in the phenomena from being studied in 
a single work. Any attempt to construct a 'first principles' 
model is hindered by uncertain intermediate steps relating 
the physics of the central engine to the properties of the rel- 
ativistic jet and the gamma-ray emission mechanism. Nev- 
ertheless, in this paper we argue that the magnetar model 
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is uniquely predictive. This allows us to construct a self- 
consistent model which can in principle be compared di- 
rectly with observations. Although we focus on magnetars 
formed via the core collapse of massive stars, we also apply 
our results to AIC r ^6.7[) . Our primary conclusion is that a 
remarkable fraction of GRB properties find natural expla- 
nations within the proto-magnetar model. 

1.1 Black Hole vs. Magnetar 

In the original collapsar model, IWooslevI l| 19931 ) envisioned a 
'failed supernova,' in which the energy released by core col- 
lapse is insufficient to unbind the majority of the star, such 
that a black hole necessarily forms. If the collapsing enve- 
lope has sufficient angular momentum, it accretes through a 
centrifugally-supported disk. Energy released by accretion, 
or vi a the accretion-mediated ex traction of the black hole's 
spin l|Blandford fe Znaieldll977r i, then powers a relativistic 
jet, which burrows through the star and ultimately pow- 
ers the GRB at larger radii (iMacFadven fe Wooslev 1 19991 ; 
iProga et al.ll2003l ; lMatzneij[2003l ; lMorsonv et al.ll2007h . 

The discovery that long GRBs are accompanied by 
hyper-energetic (~ 10 52 erg) SNe propelled the collapsar 
model to the theoretical forefront. However, it also proved, 
somewhat ironically, that G RB-SNe are far from the com- 
plete 'failures' envisioned by 1 Wooslevl (1993). Indeed, if the 
collapsar scenario is correct, then either (1) the BH forms 
promptly following stellar collapse and the explosion mech- 
anism associated with GRB-SNe is fundamentally different 
than that associated with the death of normal (slower rotat- 
ing) stars, which are instead powered by NS formation; or 
(2) a BH forms only after several seconds delay, due to the 
'fall-back' of material that remains gravit ationally bound 
despite a successful and energ etic SN (e.g.|Chevalierill993l: 

2oid T 



Fryer 1999; Zh ang et al. 2008; M oriva et al 



Modern core collapse simulations find that the shock 
produced at core bounce initially stalls due to neutrino 
and photo-dissociation losses (e.g. Rampp fe Jankal |2000| ; 



am pp fe Jan 

iLiebendorfer et alJl200j ; lThompson et alj|2003h ~ It has long 
been thought that neutrino heating from the proto-NS 
may revive the shock , resulting in a successful explosion 
|Bethe fe Wilsonlll985h . Recent simulations suggest that the 
neut rino mechanism m ay work for low mass progenitors 
(e.g. IScheck et al1l2006h . but higher mass stars appear more 
difficult to explode. Although multi-dimensional effects not 
captured by pr esent simulations may be a crucial missing 
ingredient (e.g. iNordhaus et al.ll2010f ). neutrinos alone may 
well prove incapable of powering ~ 10 52 erg explosions. 

GRB progenitors are, however, far from typical. Essen- 
tially all central engine models require rapid rot ation and a 
stron g, large-scale mag netic field (> 10 15 G; e.g. iMcKinnevI 
2006). These ingredients may go hand-in-hand in core col- 
lapse because differential rotation provides a source of free 
energy to power field grow th, via e.g. an a — Q. dyna mo 
in the convective proto-NS l|Duncan fe Thompsonlll992h or 
the magneto-rotational i nstab ility (MRI; e.g. Akivama et al.l 
120031 ; iThompson et all [2005). The crucial question then 
arises: Do SNe indeed fail and lead to BH formation if the 
progenitor core is rapidly rotating? or stated more directly: 
Are the requisite initial conditions for the collapsar model 
self-consistent? 

An additional energy reservoir (rotation) and means for 



extracting it (magnetic fields) make magneto-rotational ef- 
fects a more pro mising way to produce hypernovae than neu- 
trinos alone (e.g. iLeBlanc fe Wilsonll 19701 ; ISvmbalistvlll984l ; 
lArdelian et all 120051 '). Only recently, however, have simula- 
tions begun to captu re the combined effe cts of MHD and 
neu trino heating ( e.g. Burrows et al.|[2007h . 

iDessart et al .1 j2008), hereafter D08, calculate the col- 
lapse of a rotating 35M f^ ZAMS collapsar progenitor of 
IWooslev fe Hegen (2006), which they endow with a pre- 
collapse magnetic field that results in a ~ 10 15 G field 
strength when compressed to NS densities. This reproduces 
the field strength, if not the field topology, expected from 
the saturated state of the MRI. Soon after core bounce, a 
bipolar MHD-powered outflow develops from the proto-NS. 
Although the explosion is not initially successful over all 
solid angles, matter continues to accrete through an equa- 
torial disk. By accreting angular momentum, the NS re- 
mains rapidly spinning, which in turn enhances the mass 
loss from higher latitudes due to magneto-centrifugal sling- 
ing (e.g. IThompson et alj 12004 iMetzger etail 120071 ; see 
eq. |A12| ). Importantly, in the strongly magnetized model 
of D08, the wind mass loss rate eventually exceeds the ac- 
cretion rate, such that for t > 300 ms the NS mass be- 
gins decreasing. Although D08 cannot address the possi- 
bility of later fall-back, and a different progenitor angular 
momentum profile could change the conclusion, their result 
is nonetheless suggestive: a core self-consistently endowed 
with the properties required to produce a GRB may not 
leave a BH at all. The results of D08 highlight the fact that 
BH versus NS formation may not be a function of progeni- 
tor mass and metallicity alone. Delineating this dichotomy 
more definitively will, however, require addressing challeng- 
ing theoretical issues, such as the precise me chanism re spon- 
sible for amplifying the magnetic field (see ISpruitll2008l for 
a discussion). 

Figure Q] is a schematic diagram of the possible effects 
of rapid rotation and strong magnetic fields on the regimes 
of NS versus BH formation as a function of main-sequence 
stellar mass M* and the initial NS rotation period Po- The 
collapse of slowly rotating, low mass stars may result in a 
normal SN with kinetic energy ~ 10 51 ergs powered by neu- 
trinos. For higher mass stars, however, neutrino-powered ex- 
plosions are less likely (or are accompanied by significant 
'fall-back' accretion) due to more massive, compact iron 
cores and higher envelope binding energies i?bmd- For these 
reasons it has been argued that stars with M+ > 25Mq leave 
BH remnants at the sub-solar met allicities that appear to 
characterize GRB progenitor s (e.g. lFrver|[l999l ; iHeger et al.l 
120031 ; lO'Connor fe OttfeoiOh . 

Above the dashed line in Figure [TJ however, the ro- 
tational energy E TOt of the proto-NS (eq. [T]) exceeds the 
binding energy of the stellar envelope, where 



(i/2)m 2 



., 1n 52 I M ns 

3 x 10 ergs 



I AM, 







12 W \ms) ' [ ' 



and I = (2/5)M ns _R 2 s , M na , i? ns , and Q = 2n/P are the 
NS moment of inertia, mass, radius, and rotation rate, re- 
spectively. We have defined -Ebin d exterior to 1.8Mq, as cal- 
culated bvlDessart et all (|2010I ) from the stellar profiles of 
I Wooslev et all l|2002l ). Although the efficiency with which 
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E t ot couples to the SN shock depends on uncertain details 
during the first few hundred milliseconds after core bounce, 
if E TO t > -Ebind then a NS remnant could in principle result, 
even for very massive stars. The hypothetical boundary be- 
tween NS and BH formation based on the above discussion 
is shown with a solid line in Figure [T] We note that there 
is indeed evidence that some Galactic m agnetars may hav e 
stellar progenitors with masses J> 40Mq (|Muno et al.ll 2006), 
although (consistent with Fig. P this do es not exclusively 
appear to be the case (|Davies et al.ll2009h . 

If an MHD-powered SN does not leave a BH, then 
a rapidly spinning, strongly magnetized NS (a 'proto- 
magnetar') may instead remain behind the outgoing SN 
shock. The rotational energy E TOt > 10 52 ergs of a magnetar 
with Po ~ 1 ms is more than sufficient to power most long 
GRBs. However, not all of this energy is available to pro- 
duce high energy emission; a fraction of E Iot , for instance, 
is expended as the jet emerges from the star or is used to 
power an accompanying hypernova (dashed line; Fig.[T}. The 
right axis in Figure [T] shows the magnetic field strength B cq 
that would be generated if the magnetic energy in the dipole 
field is ~ 0.1% of E la t (eq. [3]). A dot-dashed line shows the 
minimum rotation rate required to produce a classical GRB 
from a magnetar with a field strength Pdi P , based on the 
model presented in The conditions for a hypernova and 
a GRB from a proto-magnetar are thus remarkably similar. 



description of the individual model stages described above 
using an extensive parameter study of wind models. 

Although the site and mechanism of prompt GRB emis- 
sion remain uncertain, in §3] we calculate the light curves 
and spectra within two emission models. Depending on the 
means and efficacy of the jet's acceleration ( H4.ip , GRB 
emission may be powered by the dissipation of the jet's 
Poynting flux directly ('magnetic dissipation'; i]4.3[) near or 
above the photosphere; and/or via 'internal shocks' within 
the jet at larger radij^ f ij4.4p . As Figure [2] makes clear, self- 
interaction in the jet is inevitable because ao — and hence the 
jet speed— increase monotonically as the proto-NS cools. 

After t ~ 30 — 100 seconds, ao increases even more 
rapidly as the proto-NS becomes transparent to neutrino 
emission. Because magnetic dissipation and jet acceleration 
become ineffective when ao is very large, this abrupt transi- 
tion likely ends the prompt GRB. In i|5]we address the pos- 
sibility that residual rotational or magnetic energy may con- 
tinue to power late time flaring or afterglow emission, such 
as the X-ray plateau. In Sj6]we discuss the implications of our 
results for the diversity of GRB-related phenomena, includ- 
ing very luminous GRBs ( N6.2[l . low luminosity GRBs f £|6.3f> . 
thermal-rich GRBs/X-ray Flashes ( i]6.4|l . Galactic magne- 
tars ( H6.6p , very luminous supernova f tj6.5p . and magnetar 
formation via AIC ( i]6.7jl . We summarize our conclusions in 

m 



1.2 Summary of the Magnetar Model and This 
Paper 

In this section we summarize the organization of the paper 
and orient the reader with a brief description of the model 
timeline (more details and references are provided in subse- 
quent sections). 

In §2] we present calculations of the time-dependent 
properties of proto-magnetar winds and quantify the stages 
of the proto-magnetar model. The basic picture is summa- 
rized by Figure [2] which shows the wind power E and mag- 
netization (To (maximum Lorentz factor) as a function of 
time following core bounce, calculated for a proto-magnetar 
with a surface dipole magnetic field strength Pdi P = 2 x 10 
G, initial spin period Po = 1.5 ms, and magnetic obliquity 
X = 7r/2. Changes in the wind properties with time are 
driven largely by the increase in oo(t) as the proto-NS cools. 

Within the first few hundred milliseconds following core 
bounce, a successful SN shock is launched by neutrino heat- 
ing or MHD forces (Stage I). Soon after, a wind heated by 
neutrinos expands freely from the NS surface into the cav- 
ity evacuated by the outgoing shock. The wind is initially 
non-relativistic (<to < 1) because the neutrino-driven mass 
loss rate is high (Stage II). However, as the proto-NS cools, 
ao increases to ;> 1 and the wind becomes relativistic (Stage 
III). The wind is collimated by its interaction with the star 
into a bipolar jet, which breaches the stellar surface after 
~ 10 seconds. After jet break-out, the relativistic magne- 
tar wind is directed through a relatively clear channel out 
of the star and the GRB commences (Stage IV; Aver- 
aging over variability imposed by e.g. interaction with the 
jet walls ( H4.2p , the time evolution of the power and mass- 
loading of the jet match those set by the magnetar wind at 
much smaller radii. In Sj3] we provide a more quantitative 



2 PROTO-MAGNETAR WINDS 

In this section we present calculations of th e time-dependent 
properties of magnetized proto-NS winds l|Thompson et al.l 
l2004l;lMetzger et alj|2007h . In ffTTl we summarize the model , 
which is similar to that presented in iMetzger et ail (|2007T ) 
but includes additional details not addressed in previous 
work. Our results are presented in £12.21 



2.1 Evolutionary Wind Model 

2.1.1 Model Description 

The two most important properties of the proto-magnetar 
wind are the mass loss rate M and the energy loss rate, 
or wind power, E. The wind power contains kinetic and 
magnetic (Poynting flux) components: E — Pkin + Pmag- A 



related quantity, determined from M and P m a g , 
magnetization^ 



Mc 3 



is the wind 



(2) 



1 In this paper we define 'internal shocks' as those resulting from 
the interaction between the magnetar jet and the accumulated 
(slower) shell of material released at earlier tim es. This is in con- 
trast to the standard internal shock model (e.g. |Rccs & M eszarosl 

which invokes the singular interaction between shells with 
similar properties released immediately after one another. As we 
discuss in i|4,4l the former dominate the latter in the magnetar 
model because the mean Lorentz factor of the jet increases mono- 
tonically in time. 

2 Note that this definition may differ from that used else- 
where in the literature. In particular, what we define as ao 
is s ometimes referred to as the 'baryon loading' parameter 
(e.g. iDrenkhahn fc Spruit|[2002l) . 
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Time (s) 

Figure 2. Wind power E {right axis) and magnetization crp {left axis; eq. [2]) of the proto-magnetar wind as a function of time since 
core bounce, calculated for a neutron star with mass M nB = 1.4Mq, initial spin period Pp = 1.5 ms, surface dipole field strength 
^dip = 2 X 10 15 G, and magnetic obliquity % = t/2. Stages denoted I.— V. are described in detail in [[3] 



where 0. is the NS rotation rate, (j> = B r r 2 is the magnetic 
flux t hreading the open magnetosphere divided by 4ir stera- 
dians (Michel 1969), and B T ~ the poloidal field strength. As 
shown in Appendix [XI <f> is directly related to the Poynting 
flux _E m ag (eqs. [A1| . |A3] ). The magnetization is important 
because it delineates non-relativistic (<7q < 1) from relativis- 
tic (crp > 1) outflows and affects the asymptotic partition 
between kinetic and magnetic energy in the wind. In partic- 
ular, in relativistic outflows most of the wind power resides 
in Poynting flux (£ mag £"kin) at the fast magnetosonic 
surface. The value of &o in this case crucially affects the ef- 
ficiency with which the jet may accelerate and dissipate its 
energy f i|4.1[) and is approximately equal to the outflow's 
maximum achievable Lorentz factor r max « E/Mc 2 ~ ao- 

In Appendix [X] we describe in detail how E, M, and 
ao are determined in magnetized proto-NS winds. To briefly 
summarize, mass loss during the first t ~ 30 — 100 seconds is 
caused by neutrino heating in the proto-NS atmosphere. As 

5/3 10/3 

a result, M oc LJ t v depends sensitively on the neutrino 
luminosity L v and the mean neutrino energy e v during the 
Kelvin-Helmholtz cooling phase (eq. [A8| ~). In most cases we 
t ake L v {t) and e v (t ) from the proto-NS cooling calculations 
of I Pons et alj l| 19991 ) (see Fig. lAl|) . but modified by a 'stretch 



factor' 7/ s (defined in eq. [All] ) that qualitatively accounts 
for the effects of rotation on the cooling evolution. 

We assume that mass loss from the proto-NS occurs 
only from portions of the surface threaded by the open mag- 
netic flux. We assume a dipolar magnetosphere, bounded 
by the bundle of 'last-closed' field lines which intersect the 
'Y' point radius in the magnetic equator (Figure [3] is an 
illustration of the relevant geometry). We determine the 
dependence of the Y-point radius on the wind properties 
using results from the axisymmetric MHD simulations of 
iBucciantini et al] (|2006l ) . which span the up < 1 to gp > 
1 tran sition. Using numerical results from iMetzger et al.l 
(2008), we further account for the enhancement in M that 
occurs due to magneto-centrifugal forces in the heating re- 
gion. This effect is most important when the NS is rotating 
very rapidly (P < 2 ms) and the magnetic obliquity is large, 
such that the polar cap samples regions near the rotational 
equator. After t = t„_thm ~ 30 — 100 seconds, the proto-NS 
becomes transparent to neutrinos, which causes L v and e„ to 
decrease sharply (Fig. lAl| ). Once neutrino heating decreases 
sufficiently, other processes (e.g. 7 — B or 7 — 7 pair produc- 
tion) likely take over as the dominant source o f mass- loading 
iHibschman fc Aronsll200ll ; lThompsorj|l200Sh and the wind 
composition may change from baryon- to pair-dominated. 
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(3) 



Figure 3. Geometry of magnetized proto-neutron star winds. 
The neutron star radius R nB is initially large (> 20 km) following 
the launch of the supernova shock, but decreases to its final value 
R nB 12 km in a few seconds (Fig. lAlll , The neutron star rotates 
at an angular velocity SI = 2ir/P about the vertical axis, where P 
is the rotational period; the light cylinder radius is i?L = C /Sl — 
50(P/ms) km. The magnetic dipole moment = Bdip^ns makes 
an angle \ with respect to the rotation axis. The angle O pen 
defines the size of the open magnetosphere on the neutron star 
surface. The magnetosphere is closed at angles > O pen/2 from 
the magnetic pole, while field lines with 9 < 9 opcn /2 form an 
'open' or 'wind' zone along which matter may escape to infinity. 
The size of the open zone affects both the spin-down rate and the 
mass loss rate from magnetized proto-neutron star winds. The 
bundle of last closed field lines intersects the magnetic equator 
at the 'Y' point radius Ry- Ultra-relativistic, force- free winds 
(cto 3" 1) have Ry ~ i?L> while less magnetized winds in general 
have R Y < -Rl (see flA2l and Fig.HJ. 



Lacking a predictive m odel for M at late t i mes, we assume 
that M scales with the lGoldreich fc Julianl (|l969h flux for a 
fixed value of the pair multiplicity = 10 6 . Our conclu- 
sions are fortunately insensitive to this choice (see Sj5|. The 
full expression for M is given in equation (|A15J> . 

Proto-magnetar winds are magnetically-driven 
throughout most of their evolution. When the wind is 
non-relativistic, its speed at the fast surface is v^, ~ Oq c, 
the wind power is E oc a^ 3 M oc M 1//3 and B mag = 2E^ ln 
l|Lamers fc Cass inclli 1999). For relativistic winds E oc o-qM 
is approximately independent of M, and E ma , & S> -Ekm at 
the fast point. Indeed, in the limit that o~o S> 1 we assume 
that E approaches the force- free spin-down rate (jSpitkovskvl 
2006), which depends only on <fi and SI. Even for relatively 
large (but finite) values of crrj, however, spin-down occurs 
more rapidly than in the force-free case because the 'Y' 
point radius Ry resides inside the light cylinder (see Fig. [3]). 
The full expression for E is given in equation (|A5[) . 



2.1.2 Spin-Down Evolution and Initial Conditions 

Proto-magnetar winds are magneto-rotationally powered 
throughout most of their evolution. The NS thus loses angu- 
lar momentum J = ISl to the wind at the rate J = —E/Sl. 
Neglecting mass loss (a good approximation), the rotation 



rate SI evolves according to 

SI _ 2R nB 2E 
SI R ns E TO t 

where E Iot is the NS rotational energy (eq. [I]). In equa- 
tion (0 we neglect angular momentum losses due to grav- 
itational waves, which become important if the NS is suf- 
ficiently aspher ically distorted by its st r ong interior mag- 
netic field (e.g. ICutlerll2002l ; lAronsll2003l ; IStella et al]|2005l ; 
IPairOsso et al.l l2009l ). This is a good approximation pro- 
vided that either the magnetic obliquity is small or the 
interior magnetic field is less than ~ 100 times stronger 
than the outer dipole field. We also neglect gravitational 
wave emission due t o non-axisymmet ric waves or instabil- 
ities (e.g. r-modes; lAnderssonl Il99ct) . although these are 
implicitly taken into account through the maximum initial 
NS rotation rate that we consider (see below). We also ne- 
glect the possib ility of late-time acc r etion onto the proto- 
magnetar (e.g. iMetzeer et all l2008l : IZhang fc Dal l2009h . 
which could affect the spin-down evolution both through 
accretion torques and by altering the geometry of the mag- 
netosphere. 

Given E and M as a function of SI and time, we solve 
equation (J3j) to obtain Sl(t), M(t), E(t), and a (t). A wind 
solution is thus fully specified by just four parameters: the 
NS mass M nB ; the 'initial' angular rotation rate Slo = 2n/Po; 
the surface dipole magnetic field strength Bdi P ; and the in- 
clination angle x ('obliquity') between the magnetic and ro- 
tational axes (see Fig. [3j. Since the proto-NS is still con- 
tracting for several seconds following core bounce, fio and 
Bdip are more precisely defined as the maximum values 
that would be achieved were the NS to contract at con- 
stant angular momentum J oc R^ s M nB Sl and magnetic flujjfl 
$ oc -Bdip-Rns, respectively. 

If the magnetic field is amplified on a timescale com- 
parable to the duration of the NS cooling epoch (e.g. via 
linear field winding), the assumption of a fixed dipole flux 
may be a poor approximation. On the other hand, if field 
rowth occurs more rapidly v ia a convection-driven dynamo 
..Duncan fc Thompson! I1992T) or the dynamical-tim escale 
MRI (e.g. lAkivama et alj|2003l ; iThompson et al.ll2005t ), then 
the field is pr obably established - and fin ds a MHD stable 
configuratio n llBraithwaite fc Spruitl l2006n - in less than a 
few seconds (|Spruitll2008l ). In this case the assumption that 
$ is fixed may be reasonable. 

Given the uncertainty in the origin of magnetar fields, 
in general we allow both Po and Bdip to vary independently 
within their respective physical ranges (Po > 1 m s, Bdi P < 
3 x 10 16 G; see below). However, if the magnetic field is in 
fact generated from the free energy available in differential 
rotation, then a relationship between Bdip and Po of the 
form 



1/2 ( R ™ Y 1/2 

U2km/ \ms) 



(4) 



5diP = 10 16 G( I ^ 

could result, where we have assumed that the magnetic 
energy in the dipole field (oc B^ ip R^ s ) is a fraction eb 
of the rotational energy E TOt oc Pns-Po" 2 ( e T- H) and 

3 Note the distinction between the conserved dipole flux through 
the stellar interior 3> defined here and the open flux through the 
magnetosphere <f> (eq. [2]), which evolves in time. 
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Figure 4. Time evolution of the light cylinder radius Rl (solid 
line), Alfven radius _Ra (dot-dashed line; eq. |A2| ) , 'Y' point ra- 
dius Ry (dashed line), sonic radius R B (double dot-dashed line), 
and neutron star radius R ns (see Fig. IA1JI for the solution shown 
in Figure [2] 
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Figure 5. Same as Figure [2] but calculated for different proto- 
magnctar properties. The first two models are for Po = 1 ms 
and B dip = 10 16 G, and assume values of the magnetic obliquity 
X = tt/2 (dotted line) and X = (solid line), respectively. The 
dashed line shows a lower spin-down case, calculated for Po = 2 
ms, P dip = 10 15 G, and x = n/2. 



individually in the next section. At late times ao plateaus 
and then begins decreasing once the wind mass loss rate 
reaches its minimum value proportional to the Goldreich- 
Julian flux (eq. |A14| h Once cro > 1 force- free spin-down 
obtains, such that E asymptotes at late times to the stan- 
darcQ force- free decay E oc t~ 2 . 

Figure[S]shows three additional wind models, calculated 
for different values of Bdip, Po, and %■ The models shown 
with solid and dotted lines correspond, respectively, to high 
spin-down cases with Bdip = 10 16 G, Po = 1 ms, calcu- 
lated for different values of the magnetic obliquity \ — 
and 7r/2. The third model shown with a dashed line is a 
lower spin-down case with Bdip = 10 15 G, Po = 2 ms, and 
X = n/2. Although the evolution of E(i) and ao(t) are qual- 
itatively similar to the fiducial model in Figure [2] differences 
are apparent. Note that the higher(lower) spin-down mod- 
els achieve larger(smaller) values of E and cro, except at late 
times. Also note that at fixed Bdi p and Po, ao is larger for the 
aligned rotator (x ~ 0) than in the oblique case (x ~ t/2) 
due to the enhanced mass loss in the latter case caused by 
centrifugal 'slinging' (see eq. [A 12] and surrounding discus- 
sion). 

Table Q] summarizes the results of several additional cal- 
Po = 1.5 ms, Pdip = 2 x 10 15 G, and x — tt/2. Figure 0] culations, which explore the sensitivity of our results to vari- 
shows the time evolution of several critical radii associated ations in the proto-magnetar properties and in the adopted 
with this wind solution. NS cooling model. Our primary conclusion is that key ob- 

During the first few seconds, E rises because Q and B dip servables are most sensitive to the dipole field B dip , rota- 
increase by angular momentum and magnetic flux conserva- tion rate p o, and obliquity X - Plausible variations in the NS 
tion, respectively, as the proto-NS contracts to its final ra- mass M »> stretch parameter r] s , and the cooling model, on 
dius. On longer timescales, E reaches a maximum and then the other hand, generally result in at most order unity dif- 
decreases once the NS begins to spin down and the open ferences. For this reason we fix M ns = 1.4M and rj s = 3 
magnetosphere shrinks. The latter results because both the in the sections to follow and confine our analysis to the 3D 
spin-down and the larger wind magnetization cause R Y to parameter space (B dlp , P , x)- 
increase (see Figs. [3] and [4|. Figure [2] also shows that ao in- 
creases rapidly for the first ~ 100 seconds as the NS cools 

and the neutrino-driven mass loss rate decreases. This re- 4 Note) howeverj that the measured braking indices of Galac- 

sults in several distinct stages in the wind evolution, which tic pulsars generally diffe r from the force-free prediction 
we denote by Roman numerals in Figure[2]and are discussed fe.e. iLivingstone et al.ll2007l ; see 33}. 



that the energy in differential rotation scales with E rot . 
In our models we require that Po > 1 ms because this 
is th e allowed range of stable proto-NS rotational periods 
(e.g. IStrobel et alJll999l ). This maximum rotation rate may 
be enforced in practice by the efficient loss of angular mo- 
mentum incurred by very rapidly spinning NSs to MRI- 
generated turbulenc e or waves radiated by nonaxisymmet- 
ric instabilities (e.g. 



IWheeler k, Akivama 



5 or waves radiated Dy nonaxisymmet- 
Thompson et al.ll2005l ; lOtt et alJ kOQS: 



20071 ) . We furthermore only consider 



models with Pdip < 3 x 10 16 G because although fields up 
to ~ 3 x 10 17 G are in principle possible if eb ~ 1, sta- 
ble magnetic configurations generally require a total field 
strength whic h is larger than the dipole compo nent by a fac- 
tor ^ 10 (e.g. lTavlerlll973l ; iBraithwaitel 120091 '). In addition, 
our assumption that the magnetic field does not affect the 
neutrino-driven mass loss rate is invalid for Pdi D > 3 x 10 16 
G (see Appendix [Aj. 



2.2 Results 

The results of our calculations are summarized in Figures 
[21 -T51 and Table [T] As already discussed, Figure [5] shows 
the wind magnetization ao(t) and power E(t) as a function 
of time since core bounce, calculated for M ns = IAMq, 
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3 STAGES OF THE PROTO-MAGNETAR 
MODEL 

In this section we describe the stages of proto-magnetar 
wind evolution and quantify their relationship to GRB phe- 
nomenology. Our discussion is guided closely by Figures 

I. Pre-Supernova/Thermally-Driven Wind 

Oo < 10" 3 ; t < fewxlOO ms) 

Simulations of core collapse fail to produce a prompt 
explosion, suggesting that the proto-NS continues to ac- 
crete for several hundred milliseconds before a delayed ex- 
plosion occurs. The proto-NS forms hot and its initial ra- 
dius exceeds ~ 30 km. Since magnetic forces are unlikely to 
be dynamically important yet, an explosion at t his stage 
would be neutrino-driven (|Bethe fc Wilson! Il985l ). If this 
'standard' scenario applies, thermal pressure is initially re- 
sponsible for accelerating the neutrino-heate d wind into the 
cavity behind the outgoing SN shock (e.g . Burrows et alj 
1 19951 : iQian fc Wooslevlll996i ; iRoberts et alJl201Ch . 

However, as already discussed, hypernovae are proba- 
bly not powered by neutrinos alone. For proto-magnetars 
the field is eventually amplified to a dynamically-relevant 
strength. If this field mediates the transfer of a signifi- 
cant fraction of the rotational energy (> 10 52 ergs) to the 
SN shock, the resulting explosion would indeed be hyper- 
cncrg etk (gmjlEl For an MHD-powered SN, the neutrino 
wind is thus magnetically-driven from its onset. The division 
between thermally- and magnetically-driven winds occurs 
at a critical magnetization ao ~ 10 -3 , because above this 
value the asymptotic speed of a magnetically-driven wind 



exceeds the speed v 



0.1 c obtained via 



1/3 

^oo — lTq C 

thermal acceleration alone (jMetzger et al] |2007), 

II. Magnetically-Driven, Non-Relativistic Wind 

(10~ 3 < (To < 1; fewxlOO ms < t < few s) 

Regardless of whether the SN itself is powered by 
thermal or magnetic forces, the neutrino wind becomes 
magnetically-driven (ao > 10~ 3 ) less than a second later. 
Because the neutrino luminosity L v is still large at these 
early times (Fig. IA1[I . the wind mass loss rate M remains 
high. Though powerful at this stage, the outflow is thus still 
non-relativistic (oo <C 1). Non-relativistic magnetiz ed winds 
are e fficiently self-collimated by hoop stresses (e.g. ISakurail 
Il985h . The proto-magnetar wind thus forms a bipolar jet, 
which catches up to the slower SN shock and begins boring 
a collimated cavity into the unshocked star. 



III. Magnetically-Driven, Relativistic Wind 
(Pre-Breakout) 

(1 < cr < 10 - 100; few s < t < t bo ) 




Figure 6. Contours of the wind magnetization at jet break-out 
* = *bo = 10 s, as a function of the magnetic field strength -Bdip 
and initial rotation period Pq of the magnetar. Solid and dotted 
lines show calculations assuming magnetic obliquities \ = and 
X = respectively. 



As the NS continues to cool, <tq exceeds unity 
within a few seconds and the wind becomes rela- 
tivist ic. Self-collimation fails in ultra-relativistic outflows 
(e.g. iBucciantini et all [2006). The wind power thus be- 
comes concentrated at low latitudes, where it collides with 
the slowly-expandin g SN ejecta and forms a hot 'proto- 
magnetar nebula' l|Bucciantini et aL I 120071 ). As toroidal 
flux accumulates in the nebula, magnetic forces - and 
the anisotropic thermal pressure they in duce - redirect 
the equatorial outflow towards the poles dBegelman fc Lil 
19921; iKonigl fc Granotj[2002l; [Uzdenskv fc MacFadvenll2007l; 
Bucciantini et all 120071 . I20081 . 120091 ; iKomissarov fc Barkovl 
20071) . Stellar confinement thus produces a mildly relativis- 
tic jet, which continues drilling a bipolar cavity where the 
earlier non-relativistic outflow left off. 

The jet propagates through t he star at a significant 
fract ion /3 of the sp ee d of li ght (e.g. MacFadven fc Woosleyl 



19991; lAlqv et alj |2000|; iRamirez-Ruiz et al 



2002 



Izhang et al.l 120031 ; I Morsonv et al.l 120071 ; IBucciantini et al.l 
l200ct ). such that it 'breaks out' of the stellar surface of 
radius i?* on a timescale 



t ho w R*/Pc ~ 7(^/10" cm)(/3/0.5)~ 



(5) 



Although the precise value of tbo will in general depend on 
both the properties of the jet and star, in what follows we 
assume a fixed value tbo = 10 seconds. Although this is 
a reasonable estimate for moderately powerful jets, weaker 
jets could require significantly longer to reach the surface. 
Below a critical jet power E <. J? m in ~ 10 4 8 erg s _1 , both 



hydro dynamic (e.g. [Wooslev fc ZhangH20"07r i and MHD out- 
flows (IBucciantini et all 120091) may fail to produce stable 
clean jets (e.g. lMatznerll2003l ) which may instead be 'choked' 
inside the star, resulting in little direct electromagnetic ra- 
diation (see >16.5|l . 



5 Note also that the large temperatures behind the shock pro- 
duced by such an energetic explosion will result in a large yield 
of 56 Ni. 



IV. Magnetically-Driven, Relativistic Wind (GRB) 

(10 - 100 < cr < 10 4 ; t bo <t< t cnd ) 
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After the jet breaches the stellar surface a rela- 
tivel y clean opening is so on established through the star 
(e.g. iMorsonv et all [20071 ) . Simulations suggest that after 
this point the power and mass loading of the jet reflect, in 
a time- and angle- averaged sense, the values of E(t) and 
M(t ) set by the proto-magnetar wind at much sm aller radii 
(e.g. IBucciantini et al.ll2009l ; IMorsonv et"al]|2010h . 

Figure [6] shows contours of the wind magnetization ao 
at break-out (t = tbo = 10 s), as calculated using a grid of 
wind models spanning the physical range of magnetar pa- 
rameters Bdip and Po for two values of the magnetic obliq- 
uity X = 0, 7r/2. Note that high spin-down magnetars (up- 
per left corner) produce outflows that are ultra-relativistic 
at break-out, i.e. ao\t ho ^ 10 — 100. 

Over the next tens of seconds ao increases from o"o|t bo to 
> 10 4 (Figs. [2] and [5]), resulting in ideal conditions for high 
energy emission. Assuming that the wind is collimated into a 
jet with a half-opening angle 9j , the 'isotropic' jet luminosity 
-Eiso is larger than the wind power E by a factor , where 
/ b ~ Of/2 is the beaming fracti on fR hoads 199^). Using 
axisymmetric MHD simulations, IBucciantini et all (|2009l ) 
found #j ~ 5—10° for a magnetar with Bdip ~ 3 x 10 15 G and 
Po ~ 1 nis, values consistent with the typical o pening an- 
gles inferred from GRB a fterglow modeling (e.g. I Frail et ail 
l200ll ; lBloom et alj|2003ah . 

Although the more general dependence of #j ot on the 
properties of the magnetar and stellar progenitor has not 
yet been determined, some insight is provided directly 
from observations. By combining the well-known correla- 
tion between the peak energy of the prompt emission spec- 
trum ig D eak and th e isotropic energy P iso , P pca k oc E° s * 
|Amati et alj|2003 ) with the correlation P poa k oc P 7 7 be- 
tween B Dea k and the b eaming-corrected energy P 7 = /b-Eiso 
(Ghirlanda et al. 20(3), we obtain the empirical relationship 



(cf. iNava et al. 



2006) 



/ b « 2 x 10" 



E~, 



10 51 ergs 



-3/4 



3.3° 



E~< 



10 51 ergs 



-3/8 



In what follows we assume for simplicity a fixed beaming 
fraction A, = 2 x 10 -3 , but we return to an implication of 

3/4 



the correlation /b oc in 

To produce high energy emission the jet must both ac- 
celerate to a high Lorentz factor Fj ~ ito > 1 and dissipate 
much of its bulk energy internally. Both of the emission mod- 
els that we consider in 3H magnetic dissipation and internal 
shocks, predict a characteristic emission radius where most 
dissipation occurs P 7 = P mag and B 7 = Ri s , respectively, 
that increases with time. Here P mag and Ri s are the radii at 
which magnetic dissipation peaks and internal shocks occur, 
respectively (see below). Whether photons escape the emis- 
sion region at a given epoch depends on the location of P 7 
with respec t to the radius of the Thompson photosphere of 
the jet (e.g. iGiannioalbOOot ) 



Rph 



Eiso^c 

87rc 3 <7, 



(7) 



where re ca is the Thomson opacity and we have assumed 
efficient acceleration, i.e. Tj ~ ao 3> 1 ( N4.1|) . 




o: 10 ,4 r 



100 



Time (s) 



Figure 7. Photosphere radius R p ^ (solid line; eq. [7]), inter- 
nal shock radius _Ri s (dashed line; eq. |B3) ^. and the 'satura- 
tion' radius at which magnetic dissipation peaks -R mag (dotted 
line; eq. |10j ~) in the proto-magnetar jet as a function of time 
since core bounce, calculated for the model shown in Figure [2] 
The jet breaks out of the star at the time t = t^ = 10 sec- 
onds. At times tbo < t < £thin,m a g(tthin,is) magnetic dissipation 
(internal shocks) occur below the photosphere and the resulting 
emission will be thermalized (Stage IVa). By contrast, at times 
i > *thin,m a gi tthin.is emission occurs in an optically-thin environ- 
ment and may be non-thermal (Stage IVb). The end of the GRB 



is defined as when R n 



R is (Stage V). 



10" 



(6) io 1: 



t 



'hin.mag — 10s 



■t, 



"lin.is — 10s 



thin 



™ag ~ 20 s 



"iin, is = 20 S 



,. 'hin,mag = 40 s 
thin.i, — 4ID ■> 




Figure 8. Contours of the time after core bounce tthin when the 
jet becomes optically thin to emission at the magnetic dissipation 
radius -f? mag (tthin, m ag ; dotted line) and the internal shock radius 
(tthin is! solid line) as a function of magnetic dipole field strength 
-Bdip an d initial rotation period Po, calculated for \ = tt/2. Jets 
from lower field magnetars are optically thick at break-out (i.e. 
tthin > 'bo = 10 s), potentially resulting in a short-lived phase of 
dim quasi-thermal emission (Stage IVa). By contrast, jets from 
magnetars with stronger fields (upper diagram) have t t hin < tbo 
and may dissipate their energy in an optically-thin environment 
immediately after break-out (Stage IVb), thereby skipping Stage 
IVa entirely. 



IVa. Quasi-Thermal, Photospheric Emission 
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Figure 9. Contours of the fraction of energy released in thermal 
emission / tn = P tn /(&y + -Eth) during the GRB phase, where 
i?th is defined in equation {5J , as a function of surface dipole field 
-Bdip and initial rotation rate Po for magnetic obliquities x = f/2 
(solid line) and \ = (dotted line). 

(tba < t < tthini &y < Rph) 

Figure[7]shows the time evolution of the photosphere ra- 
dius J? p h and the radii at which internal shocks (Bi s ; eq. |B3| ) 
and magnetic dissipation (B mag ; eq. [10]) occur, calculated 
for the fiducial model shown in Figure [2] Just after break- 
out, magnetic dissipation and internal shocks occur below 
the photosphere, i.e. ii 7 = {-R mag ,-Ri s } <C -Rph, such that 
high energy emission will be partially thermalized and sup- 
pressed due to adiabatic losses. At later times, jet dissipation 
occurs in an optically-thin environment (B 7 > Rph), such 
that brighter non-thermal emission^ is more likely. 

Figure [8] shows contours of the time after core bounce 
at which _R p h = Rma.g (dotted line) and B p h = Ri B (solid 
line), respectively, as a function of Bdip and Po for x = tt/2. 
Low field magnetars (lower diagram) produce jets that are 
optically thick at break-out (i.e. i t hin > tbo ~ 10 s) and 
thus experience a phase of quasi-thermal photospheric emis- 
sion, as in the fiducial model described above (Stage IVa). 
In fact, if tthln becomes comparable to the GRB duration 
itself (cf. Fig. I13[) . a thermal-rich sub-luminous GRB or X- 
ray Flash may result instead of a classical GRB f £|6.4[) . By 
contrast, jets from strongly magnetized magnetars (upper 
diagram) dissipate their energy in an optically-thin environ- 
ment just after jet break out, thereby skipping Stage IVa 
entirely. 

Figure [5] shows contours of the fraction of the en- 
ergy released in thermal emission during the GRB phase 
fth ~ Eth/(E y + -Eth)- Here E-, is the total non-thermal 
emission during the GRB (quantified in the next section) 
and Et h is the maximum therm al energy, which we estimate 
as (e.g. iMeszaros fc Reedl200Ch 

/"tthin.is 

E th = / e T E(R ph /R iB )- 2/3 dt, (8) 

J tba 

6 Throughout this paper we define 'non-thermal' emission as a 
non-black body spectrum. This does not necessarily imply that 
the radiating electrons have a non-thermal energy distribution. 
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Figure 10. Contours of the maximum non-thermal gamma-ray 
emission B 7 in ergs as a function of -Bdip an< i P) f° r X = 
(solid line) and \ = n/2 (dotted line). We calculate P 7 as 
the total energy released by the magnetar in the time interval 
max[f. DO , tthin.is] < t < t en d times a factor e r = 0.5 to account for 
the maximum radiative efficiency. Here t^ Q is the time required 
for the jet to propagate through the star, t t h_i n is is the time af- 
ter which the outflow is optically thin at the internal shock ra- 
dius (Fig. [S), and t cn d is the end of the GRB, defined as when 
flmag = Ris (see Figs. FF1 and IT3)l . 

where the factor (R p h/ Ris)~ 2 ^ 3 accounts for adiabatic losses, 
and we have (optimistically) assumed a radiative efficiency 
e r = 0.5 (eq. [B7| ). High field magnetars (upper dia- 
gram) produce little thermal emission E t h ~ because 
the jet is already optically thin at break out (i.e. tthin.is < 
th ; cf. Fig. [8}. By contrast, somewhat lower- field mag- 
netars (middle-left diagram) have fth > 0.1 (E t h_ ~ 
10 48-50 ergs), consistent with measurements or upper- 
limi ts on quasi-thermal photospheric emission from GRBs 
fe.g.lMeszaros fc Reesll2000l ; |Ramirez-Ruiz et al.ll2002l ; iRvdel 
12003 : iGuiriec et al.ll201ol v 

IVb. Main GRB Emission 

(£thin < t < tend) R~l > Rph) 

From the time t = min[tthin, tbo] until the GRB ends at 
£ = tend (which we define more precisely below), shocks or 
reconnection occur above the photosphere and non-thermal 
emission is likely. Figure [101 shows contours of the total en- 
ergy released by the magnetar wind E~, = J e T Edt inte- 
grated over the GRB duration as a function of Bdip and 
Po for x = (solid line) and x — f/2 (dotted line), assum- 
ing a radiative efficiency e r = 0.5. Note that E-, ^ ^g 5 °- 51 
ergs across the entire range of high spin-down ('GRB ca- 
pable') magnetars. These values are consistent with the 
collimation -corrected energy released by GRBs in relativistic 
ejecta fe.g. lFrail etllll200 AlBerger et al.|[2003l ; iBloom et all 
l2003al ). 

Figure [11] shows contours of the average (energy- 
weighted) magnetization of the jet, which we define as 
cr avg = J Eaodt/ J Edt integrated over the duration of 
the GRB. High spin-down magnetars (upper left diagram) 
achieve values a avg ~ F max ~ 10 2 — 10 4 which are are con- 
sistent with observational constraint s on the GRB Lorentz 
factors (i.e. F > 100 — 1000; e.g. iLithwick fc Saril 1200 ll : 
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Figure 11. Contours of the (energy- weighted) average magneti- 
zation <T avg as a function of B<iip and Po, calculated for \ = 
(solid line) and x = t/2 (dotted line). Note that in general tx avg 
is smaller in the case of an oblique rotator (x = 7r/2) because of 
the enhanced mass loss due to centrifugal 'slinging' (see eq. IA12I 
and surrounding discussion). 

IZou fc Piranl |2010| ; IZou et alj|2010h . We caution, however, 
that although (j aV g approximately equals the jets maximum 
instantaneous Lorentz factor, for internal shocks the Lorentz 
factor of the emitting material F s is generally lower than 
cravg because the faster jet interacts with slower material re- 
leased at earlier times ( i|4.4[) . In Figure [12] we show contours 
of the (energy-weighted) mean Lorentz factor F s , avg of the 
bulk shell, from behind which internal shock emission origi- 
nates. Note that in general r s , avg is a factor of a few times 
lower than o- avg . A comparison of Figure [lOl with Figures [TT1 
and [T2] reveals a positive correlation between _E 7 and the 
mean magnetization/Lorentz factor. We discuss this corre- 
lation and its implications further in H4.3I 

V. Ultra High-cro Phase (Post GRB) 

(-^m a g ^ Ris> t ^ tend). 

As o~o continues to increase, the jet becomes less and less 
effective a t accelerating and dissip ating its ordered energy 
( P~T1 e.g. iLvubarskv fc Kirkl l200lT ) . A particularly abrupt 
jump in (Jo occurs once the NS becomes transparent to neu- 
trinos at t = t„_ tn in ~ 30 — 100 s (Fig. IA1[) . after which 
ao rises to very large values J> 10 9 . This transition likely 
ends the prompt high energy emission. Although the argu- 
ment for why t cn( j ~ ti/-thin is quite general, we can be 
concrete by defining t cn d as the time after which the mag- 
netic dissipation radius -R mag (eq. [TO]) exceeds the internal 
shock radius R[ s . For t > t en d the jet magnetization at the 
shock radius exceeds the crit ical value a ~ 0.1 above which 
strong shocks are suppressed ijKennel k, Coro niti 1984). The 
association of t cnt j with t„_ t hm explains both the typical du- 
ration of long GRBs T90 ~ 10 — 100 s and accounts for 
why the prompt ~ MeV emission declines more rapidly at 
lates times > Tgo (oc t -3 ) than the j et luminosity predicted 
by most central engine models (e.g. iTagliaferri et all 120051 : 
iBarniol Duran fc Kuniarll2009l) . 

Figure [13] shows contours of the rest-frame GRB 
duration T90, defined as the time interval within 
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Figure 12. Contours of the (energy-weighted) average Lorentz 
factor r s , aV g of the bulk shell created by internal shocks, as a 
function of Pdip and Po, calculated for \ = (solid line) and 
X = tt/2 (dotted line). 
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Figure 13. Duration Tgo in the time interval max[tboj tthin is] < 
t < *cnd during which 90% of the wind energy is released, cal- 
culated as a function of Pdip and Pq for X = (solid line) and 
X = tt/2 (dotted line). Here tbo is the jet break-out time, tthin 
is the time after internal shocks occur above the photosphere 
(Fig. and t cn( j is the end of the prompt emission (i.e. when 
P mag = Pi s ; see Fig. [7). 



max[tbo, tthin.is] < t < t cn d during which 90% of the wind 
energy is released. Note that high spin-down ('GRB capa- 
ble') magnetars have Tgo ~ 40 — 50 seconds, similar to the 
average observed rest-frame duration of long GRBs. A qual- 
itatively similar, though somewhat shorter, Tgo distribution 
results if we assume that emission begins at tthin, m ag (mag- 
netic dissipation) rather than tthin, is (internal shocks). The 
true predicted (rest frame) GRB duration distribution will 
of course be broader than suggested by Figure[l3]because we 
have not taken into account variations in the timescale for 
jet break-out tbo (eq. [5]) and neutrino transparency t„_ t hm, 
the latter of which depends on the NS mass and rotation 
rate. Realistic variations in tbo, Tj s , and M ns will undoubt- 
edly broaden the rest-frame Tgo distribution by factors of a 
few as observed (see Table [1] for examples). 



© ???? RAS, MNRAS OOO.HT4291 



12 B. D. Metzger, D. Giannios, T. A. Thompson, N. Bucciantini, & E. Quataert 



Except in the case of very luminous GRBs ( Ej6.2[) . most 
of the magnetar's initial rotational energy remains when the 
prompt emission ends. Though not released as gamma-rays, 
this residual energy may be dissipated at later times or larger 
radii and hence may contribute, for instance, to the GRB X- 
ray afterglow. In iJS]we discuss emission during the late-time 
high- <7o phase. 



4 GAMMA-RAY BURST EMISSION 

In this section we calculate the emission during the prompt 
phase (ibo <J t <J £ e nd; Stage IV). We begin with a discussion 
of the mechanisms for jet acceleration ( H4.ip and variability 
( £14.21) and then present calculations of the gamma-ray emis- 
sion produced by magnetic dissipation ( £|4.3[) and internal 
shocks ( H4.4I) . 



4.1 Acceleration 

Energy carried by the relativistic wind is primarily in the 
magnetic field near the light cylinder radius _Rl ~ 10 7 
cm. Because GRBs originat e from ultra-relativistic outflows 
(e.g. iLithwick fc Sarill200ll ), this magnetic energy must be 
transferred to kinetic energy prior to the radii ~ 10 12 — 
10 16 cm at which the high energy emission occurs. Uncon- 
fined, time-stationary Poynting-flux do minated outflows do 
not accelerate efficiently in ideal MHD (iGoldreich fc Julianl 
Il97d ; iBeskin et alll998i ; iBogovalov fc Tsinganosll 19991 ). The 
Lorentz factor reached when accel eration slows near the fast 
magnetosonic surface Too ~ a^ 3 l|Goldreich fc Julianlll970l ) 
is much less than the maximum possible value F max ~ 
(To. Full acceleration to Too ~ F max therefore appears to 
require a combination of a differentially-collimated (non- 
monopolar ) geometry, time v ariability, or violations of ideal 
MHD (see lKomissarovll2010l for a recent review). 

At small radii the wind is concentrated in the ro- 
tational equator. On larger scales the outflow is redi- 
recte d into a bipolar jet by i ts interaction with th e 
star (|Komissarov fc Barkovl 120071: iBucciantini et aTH2 009). 
Anal ytic (e.g. IVlahakis fc Konigll l200ll; I Nara van et al.l 
20071) and numer ical llKomissarov et alJ 12001 I2OO9F 



Tchekhovskov et al ] 120091 . l20ld ) calculations show that if 



the jet is confined into a parabolic shape, additional acceler- 
ation is possible due to 'equilibrium collimation'. However, 
the maximum Lorentz factor that can be achieved in this 
manner is Too ~ l/9j ~ 10 (eq. [6]) because only while 
r#j <J 1 does the jet remain in lateral causal contact. Al- 
though an additional boost of acceleration (by a factor <C 
10) may occur as the jet emerges from the stellar surface 
ijKomissarov et al.l [20091 : iTchekhovskov et aD 120091 ). reach- 
ing Too J> 10 2 and simultaneously achieving high conversion 
efficiency of magnetic to kinetic energy appears difficult via 
collimation alone. 

A time-dependent flow can also produce acceleration. 
In the so-called 'as trophysical p lasma gun' or 'magnetic 
rocket' mechanism (Contopoulos 1995; iGranot et alJl201ol : 
iLvutikov fc Listeril2010l : lLvutikovll2010l ). a high-00 magnetic 
pulse of finite width ex pands into a low er density medium 
('vacuum'; see, however. lLevinsonll2010l ). As the shell prop- 
agates, it 'self-accelerates' via magnetic pressure gradients 
which develop as a rarefaction wave passes through the shell. 



In this case F increases tx r 1/3 (e.g. IGranot et al.ll201oh . 
similar to the magnetic dissipation model described below. 
Faster acceleration F oc r is pos sible in standard (high en- 
tropy) GRB fireball models (e.g. lGoodman| [T986) . but it re- 
mains unclear how the necessary thermalization would oc- 
cur inside the starQ especially considering that reconnec- 
tion may be sl ow in the collisional environme nt close to the 
central engine (Mc Kinnev fc Uzden skv 2010). Note that no 
ideal MHD model for jet acceleration accounts for the dis- 
sipation of energy responsible for powering the GRB, which 
must instead occur at larger radii after acceleration is com- 
plete. 

An alternative possibility for jet acceleration is 
magnetic dissipation, i.e. a break-down of ideal MHD 
( Spru it et al ] l200ll : iDrenkhahnl l2002i : iDrenkhahn fc Spruitl 
2002). One way this can occur is if the rotation and magnetic 
axes of the NS are misaligned (x > 0), such that the outflow 
develops an alt ernating or 'striped' magnetic field geometry 
ljCoronitilll990l ) on the scale of the light cylinder radius. If 
this non-axisymmetric pattern is preserved when the flow is 
redirected along the polar jet, the resulting geometry is con- 
ducive to magnetic reconnection. Magnetic dissipation oc- 
curs gradually from small radii up to the 'saturation' radius 
fl mag , beyond which reconnection is complete and the flow 
achieves its terminal Lorentz factor. During this process, ap- 
proximately half the Poynting flux is directly converted into 
kinetic energy (producing acceleration) and the other half 
is deposited into the intern al (thermal) energy of the flow 
(|Drenkhahn fc Spruitll2002l ). Acceleration and emission thus 
bot h result from the s ame physical mechanism. 

IDrenkhahnl (2002) shows that the Lorentz factor of the 
jet as function of radius is given by 



<To(r/7? n 



\i/3 



it). 



r < R n 
r>R„ 



where the saturation radius is 



Cm (S0 (n^)( 



0.01 



(9) 



(10) 



and e <; 1 parametrizes the reconnection speed v T — ev\, 
where va — c is the Alfven speed. In our calculations we 
assume that e = 0.01, independent of radius or jet proper- 
ties. This value is motivated by recent work finding a re- 
connection rate of this order due t o secondary tearing in - 
stabilities in the current sheets (e.g. lUzdenskv et al.ll2Q10l ). 
even in highly collisional environments, that characterizes 
the jet close to the central engine. On the other hand, at 
larger radii (yet still well below the nominal saturation ra- 
dius), reconnection may occur in the collisionless r egime, 
such that faster reconnection is al so likely (see e.g. lAronsi 
|200S| . iMcKinnev fc Uzdensky||2010l for specific physical dis- 
sipation mechanisms). 



4.2 Variability 

Although GRBs ar e variable on timescales do w n to fractions 
of a millisecond ijSchaefer fc Walker! 1 19991 : IWalker et al.1 



7 One possibility is if instabilities act within the jet to randomize 
the magnetic field, such that it behaves as a 7 = 4/3 relativistic 
gas llHeinz fc Begelmanllioool; iGiannios fc Spruij|2006h . 
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Figure 14. Bolometric GRB luminosity due to magnetic dissipation (dashed line) and internal shocks (dotted line) as a function of 
observer time t b si calculated for a proto-magnetar with -Bdip = 2 X 10 15 G, Po = 1.5 ms, and \ = tt/2 (Fig. [2}. For internal shocks we 
assume that e mag = 0.5 and e e = 1 (see text for definitions). The isotropic power of the jet i?i so is shown for comparison with a solid line 
and is calculated assuming a beaming fraction /b = 2 X 10 -3 . The times when the jet becomes Thomson thin to emission from internal 
shocks and magnetic dissipation arc marked with diamonds. Although emission is suppressed at early times due to adiabatic losses, at 
times t 3> tthin the radiative efficiency of both magnetic dissipation and shocks approaches ~ 1/2 (cq. |B7| ). 



2000) , most Fourier p ower is concentrated on a c harac teristic 
timescale ~ 1 second (|Beloborodov et alJll998l . l200Cf ). GRB 
variability may be related to the emission mechanism itself, 
or it may reflect re al variations in the power and mass load- 
ing of the jet (e.g. MacFadyen fc Wooslevlll99S ; lAlov et all 
l200d : lMizuta fc Alovll2009l ; iMorsonv et alJboioh . 



There are several potential sources of variability in 
proto-magnetar outflows. Sporadic changes to the magne- 
tosphere could modulate the magnetar wind properties on 
short (< millisecon d) timescales due to rec onnection near 
the light cylinder (|Bucciantini et aTl [2006) or on longer 
timescales due t o neutrino heating in the closed zone 
|Thompsonll2003l ). Longer timescale variability could also 
be imposed on the outflow as it propagates to the stellar 
surface, due to instabilities associated wit h the termination 
shock ( s) in the proto-ma gnetar nebula (|Bucciantini et all 
120091 ; ICamus et all I2009T I or at larger di stances as the 
jet p r opaga tes through the stellar envelope l|Morsonv et al.l 
120071 . |2010| ). The latter possibility is particularly promis- 
ing because the sound crossing time across the jet near 



the stellar radius is in fa ct ~ 1 second (e.g. IMorsonv et al.l 
l20ld ; lLazzati et alJlioiol ) and might not evolve appreciably 
throughout the burst, a fact co nsistent with observations 
(R amirez-Ruiz fc Fenimordll999h . 



The time-averaged wind properties calculated in §271] 
(E, M, and 00) do not account for any of the variability dis- 
cussed above. In fact, given the stochastic nature of GRB 
emission, it seems unlikely that any model will be capable 
of predicting the detailed light curve of individual bursts. In 
our calculations below, we instead focus on predicting the 
time-averaged high energy emission over timescales of sec- 
onds or longer, which may be usefully com pared with inte- 
grate d GRB light curves and spectra (e.g. iMcBreen et al.l 
120021 ). We nevertheless emphasize that variability affects 
the observed emission differently depending on the emission 
model. Magnetic dissipation, for instance, occurs at rela- 
tively small radii, such that variability is directly encoded 
in the emitted radiation. Variability from internal shocks in- 
stead manifests indirectly through the effects of subsequent 
collisions at larger radii. 
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4.3 Emission from Magnetic Dissipation 

iDrenkhahn fe Spruitl (|2002T ) make specific predictions for the 
rate that magnetic energy is dissipated with radius. How- 
ever, reconnection can in principle energize particles in a 
variety of ways. Reconnection can lead to plasma heating 
and acceleration in localized region^f] (e.g. current layers). 
Alternatively, reconnection may drive b ulk motions in th e 
jet that excite Alfvenic turbulence (e.g. iThompsonl 1 1994 ) ■ 
which cascades to small scales and heat s larger v o lumes 
in th e plasma. We follow the model of IGianniosl (2006, 
2008), who assumes that the dissipated energy heats the 
pla sma smoothly throughout the flow (slow heating m odel; 



!Ug£ 

see lGhisellini fe Celottill 19991 : IStern fe Poutanenll2004I ) . Sim- 
ilar qualitative conclusions would, however, result from any 
model that invokes localized modest particle accelerat ion 
clos e to the p h otosp here fe.g. lLazzati fc Begelmanll2010T ). 

IGianniosl (|2008l ) shows that energy dissipated at large 
Thomson optical depths is thermalized, such that a por- 
tion emerges through the photosphere with a peak at 
MeV energies (cf. Goodman 19861; Meszaros fc Reea 



20001; iRamirez-Ruizil^OOSl ; iPe'er et al.ll2006l ; IGianniosl 120061 ; 
Beloborodovll2010l ). At times when i? mag J> -Rph most of the 
Poynting flux is dissipated near or above the photosphere 
and the equilibrium temperature of the electrons exceeds 
the radiation temperature. Inverse Compton scattering of 
the photons advected outwards with the flow then results in 
power-law emission with a flat spectral slope E ■ Le oc E 
above the thermal peak. Larger radii in the flow are heated 
to yet higher temperatures, resulting in an additional com- 
ponent of synchrotron and synchrotron-self- Compton emis- 
sion at lower frequencies (i.e. optical, UV, and X-ray bands). 
This softens the spectrum below the MeV peak close to the 
observed E ■ Le oc E 1 value. 

Figure [14] shows the bolometric (isotropic) luminosity 
due to magnetic dissipation -Lmag, calculated for the fiducial 
model shown in Figure [5] At late times t > tthin.mag ~ 20 s, 
magnetic dissipation occurs above the photosphere (i? m ag J> 
R ph ) and L mag = E iso /2 (|Drenkhahn fc Spruitl liobl ). At 
early times t <C t t hin,ma g when -R mag <, R P h, by contrast, 
Lmag is suppressed below Ei ao by an additional factor ~ 
0.4(i? m ag/-Rph) 2//3 due to adiabatic losses incurred between 
the dissipation radii and the photosphere. 

Figure [15] shows snapshots of the high energy spectrum, 
calculated at the times t = 15, 20, 25, and 30 seconds. The 
spectrum at t = 15 s corresponds to an early epoch when 
Rph 3> -Rmag and the dissipated energy is thermalized; the 
spectrum in this case is approximately Planckian with tem- 
perature T ~ 2 keV. Due to its low luminosity and X-ray 
peak, such a component of early thermal emission may be 
challenging to detect in an actual GRB (Fig. 0. At later 



times t > tthin.n 



20 s, by contrast, dissipation peaks 



near or above the photosphere. This results in a spectral 
peak at energy -E pea k ~ 10 2 keV (described below) and a 
nonthermal Comptonized tail that extends to increasingly 
higher energies as ctq rises and the outflow becomes cleaner 
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Figure 15. Spectral energy distributions E ■ Le of the magnetic 
dissipation model calculated at several times, t = 15 s (solid), 
t = 20 s (dotted), t = 25 s (dashed), and t = 30 s (dot- dashed), 
for the same model shown in Figures [2] 1141 and 1161 
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Figure 16. Peak spectral energy (or break energy) -E pca k as a 
function of observer time t t, s in magnetic dissipation (dashed 
line) and synchrotron internal shock (dotted line) models, calcu- 
lated for the model shown in Figures [2] and II II 



with time. Also note the component of synchrotron emission 
at softer X-ray/UV wavelengths, which increases in relative 
importance to the Comptonized gamma-rays as the jet mag- 
netization increases and dissipation peaks at larger radii. 

The magnetic dissipation model predicts a spectral en- 
ergy peak i? P eak (or brealjf]) similar to the observed Band 
spectrum peak ~ few hundred k eV and which is relatively 
insensitive to the jet properties. iGiannios fc Spruitl (|2007l ) 
show that to good approximation 



£peak,mag - 270 keV 



Eu 



10 52 ergss - 



\W 2 J U0V 



8 This is the approach adopted bv lLvutikov fc Blandfordl l l2003l) . 
However, because the mechanisms responsible for particle accel- 
eration in magnetic reconnection are uncertain, it is difficult to 
make concrete predictions for the resulting GRB emission in this 
case. 



9 The E ■ Le spectrum above the break may (depending on pa- 
rameters) be slowly rising. In this case the ~ MeV 'peak' is for- 
mally a break. 
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for R n 



> R 



ph • 



Figure [16] shows -E pe ak as a function of time for the fidu- 
cial model shown in Figures [14] and 1151 Note that although 
Speak rises rapidly at times t <J tthin.mag (when the luminos- 
ity is highly suppressed), E paa ^ is relatively constant during 
the GRB itself, increasing from ~ 200 keV to ~ 400 keV be- 
tween t = tthin,ma g and t = t e nd- This slow evolution results 
from the weak dependence of -E pea k on Ei ao (t) and o~o{t) in 
equation (|11[) . A rising value of S pC ak at first seems in con- 
flict with the observation that GRBs are usually inferred to 
spectrally 'soften' throughout their duration. This behaviour 
may, however, still be consistent with spectral evolution pre- 
dicted by magnetic dissipation if the synchrotron emission 
at lower frequencies begins to contaminate the soft X-ray 
bands at late times (-E P eak,ma g refers to the spectral peak of 
the Inverse- Compton emission; see Fig. I15p . 

We now consider the implications of equation ()lip for 
the population of magnetar-powered GRBs as a whole. Fig- 
ure [17] shows a scatter plot of the average magnetization 
(Tavg during the GRB (Fig. 1 1 1 p as a function of the average 
GRB luminosity L 1 = E 1 /T ao (Figs. \W\ and [13]), where we 
have included data points from all models within the range 
of magnetar parameters explored previously (1 ms < Pq < 
5 ms; 3 x 10 14 G < B dip < 3 x 10 16 G; x = and X = t/2, 
respectively). Magnetars lying on the one-parameter fam- 
ily -Bdi P oc P _1 defined by equation (Q for es = 10 -3 are 
connected with a solid line. 

Figure fTTl shows that the magnetar model predicts, with 
large scatter, a positive correlation between <T avg and L 1 . In 
particular, for the one-parameter family of solutions we find 
that o" avg oc L", where a ~ 0.5 — 1, depending on L 7 and 
X- Assuming that the GRB duration, radiative efficiency e r , 
and beaming fraction / D are similar from burst to burst, this 
correlation implies that E 1 oc (Tavg • From equation (|lll) this 
in turn leads to the prediction that i? pC ak oc E°- 2 - 0A for 
e£2_~ constant. Note that this is close to the lAmati et al.l 
(2002) relationship i5 pca k oc £*' 4 . If one furthermore drops 
the assumption that /b is constant and instead assumes 



/b oc E~ 



as motivated by the combined Amati and 



Ghirlanda relations (eq. [6]), one finds -E pea k oc E®' 3 ~ 5 , re- 
sulting in even better agreement with observations. A qual- 
itatively similar correlation is predicted between S 7 and 
the peak jet power, consistent with the related 'Yone toku' 
relation <|Yonetoku et al.l 120041 : cf. I Wei fc Gaol 120031 '). We 
emphasize that both the normalization and the slope of 
the Amati/ Yonetoku correlations are reproduced if we as- 
sume a reconnection r ate e — 10~ 2 favored by recent work 
l|Uzdenskv et al.ll2010T ). 



4.4 Emission from Internal Shocks 

If the acceleration of the jet is efficient f i]4.1[) . then a sig- 
nificant fraction of the Poynting flux is converted into ki- 
netic energy. The kinetic luminosity and Lorentz factor of 
the outflow in this case are given by Lj(t) ~ (1 — e mag )i5(i) 
and Tj(t) ~ (1 — e ma g)o"o(t), respectively, where e mag < 0.5 
is the fraction of the power radiated during the acceleration 
phase, due to magnetic dissipation f i]4.3[) . In what follows we 
assume e mag = 0.5, although e mag = would be appropri- 
ate if the magnetic energy that is dissipated is not radiated 
away or if acceleration is achieved by another mechanism. 



Because o~o ~ Tj increases monotonically during the 
GRB (Fig. [2]), slower material is released prior to faster 
material. Strong shocks will occur once the faster mate- 
rial catches up pro vided that the residual magnetization 
of th e jet is <; 0.1 (jKennel fc Coronitil [l984l ; iMimica et all 
2009). This scenario is similar to the standard internal 
shock model for GRB emission (e.g. Rees fc Meszaroslll994l ; 
iKobavashi et al.lll997l ; iDaigne fc Mochkovitchlll998f ) with a 
few key differences to be discussed below. 

Immediately after the jet breaks out of the star, the fast 
and slow ejecta have comparable energies and speeds. With 
time, however, the slow material released early accumulates 
into a common 'bulk' shell, which we characterize by its total 
rest mass M s = J t Mjdt, energy E B , mean velocity /3 B , and 

mean Lorentz factor r s = E B /M B c 2 , where Mj = Lj/TjC 2 
(see Fig. I12[) . At most times the jet's self-interaction is well 
described as a collision between the fast, variable jet and a 
slower (yet still ultra-relativistic) shell. We model this in- 
teraction using a one-dimensional kinematic model, as de- 
scribed in Appendix[B] Although this approach neglects the 
effects of pre ssure forces and the true m ulti-dimensional ge- 
ometry (e.g. IZhang fc MacFadvenllioO^ ). it provides a rea- 
sona ble first approximation to the full hydrod ynamic prob- 
lem (|Daigne fc Mochkovitclj|T99l . l2000l . T2003l ). 

Figure [2] shows the (average) bolometric luminosity 
Lis = e r Lj from shocks as a function of observer time 
tobs(l + z), calculated assuming that the fraction of the 
shock's energy imparted to electrons is e e « 1. As in the 
case of magnetic dissipation, at times t < tthm.is we suppress 
Lis by an additional factor ~ (Rph/ 'Ri s )~ 2 ^ 3 to account for 
adiabatic losses when shocks occur below the photosphere, 
where Ri B is the internal shock radius given in eq. |B3] . Note 
that although Lj decreases by a factor of ~ 6 throughout 
the burst, Lj a changes by only a factor of a few. Indeed, 
both magnetic dissipation and shock models predict that 
the bolometric luminosity should be relatively constant in 
time, a result in agreement with the appro ximately linear 
slope of cumulative GRB light curves (e.g. iMcBreen et al.l 
120021) . 

In Appendix [Bl we calculate the peak energy -E pca k,is of 
the synchrotron spectrum as a function of the jet and shell 
properties, assuming that a fraction ( e of electrons are accel- 
erated and that a fraction eb of the shock energy goes into 
generating the magnetic field (see eq. [B13| and surrounding 
discussion) . 

Figure [16] shows the evolution of £" P eak,is during the 
GRB for the fiducial model, calculated assuming e e ~ 1, 
es = 0.1, and £ e = 0.3. These microphysical parameters are 
chosen ad hoc such that i? pC ak attains a value ~ 10 2 keV 
at peak luminosity characteristic of observed GRB spectra. 
Even after this fine tuning, however, two problems remain 
for the internal shock model. First, Figure [16] shows that 
£" P eak,is increases by over three orders of magnitude dur- 
ing the burst, in contradiction with the relatively constant 
(or decreasing) peak energy measured during actual bursts. 
Although both r s and tj increase with time, Speak, is oc 
£-ip2p-4 mcreases because the jet Lorentz factor Tj in- 
creases even more rapidly (see eq. |B13] V Although a slowly- 
evolving peak energy could in principle be recovered by in- 
voking e.g. time-dependent micr ophysical parameters, fin e 
tuning appears unavoidable (e.g. IZhang fc Meszarosl f2002'). 
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Figure 17. Average magnetization <7 avg during the GRB versus the average GRB luminosity L 7 = E^/TgQ. Each point represents a 
model calculated within the range of initial spin periods 1 ms < Pq < 5 ms and surface dipole fields 3 X 10 14 G < -Bdip < 3 X 10 16 G. 
The left and right panels show calculations performed assuming the magnetic obliquity \ = and \ = tt/2, respectively. A solid line 
connects solutions lying along the one-parameter family -Bdip <x Pq 1 defined by equation assuming ejj = 10 — 3 . 



A second problem is that the variability timescale produced 
by subsequent internal collisions 8t vaT oc R Bn /2r^ oc t is pre- 
dicted to increase linearly with time, again contrary to obser- 
yations suggesting that 5t va r evolv es weakly during the burst 
(R amirez-Ruiz fc Fenimorelll999l ) F°l We conclude that syn- 
chrotron emission from internal shocks appears disfavored 
as the source of prompt emission from proto-magnetars. 



5 LATE-TIME EMISSION 

When (Jo becomes very large at late times (:> 10 5 ; Fig. [2]), 
Ruing becomes so large that, even if reconnection occurs at 
the speed of light in the co- moving frame (e ~ 1), no ef- 
ficient acceleration or dissipation occurs before the outflow 
begins to interact wit h itself or the external IS M. This is 
the 'causality limit' of iLvubarskv fc Kirkl < |200lf ). Without 
acceleration, shocks cannot occur; and without efficient re- 
connection, there can be no dissipation-powered emission. 
As we argued in fj3] this transition ends the phase of prompt 
internal emission. Similar physics occurs in the wind from 
the Crab Pulsar, for which the very high initial magnetiza- 
tion may prevent internal dissipat ion prior to the wind ter- 
mination shock at R ~ 10 cm (|Kennel fc Coronitil [l983 ; 
ILvubarskv k, K irk 2001). Emission from the nebula may be 
the result of forced reconne ction at the termination shock it- 
self lHvubarsk vibool 120051 ) or from dissipation in a striped 
wind (|Coronitilll990l ) if the pair multip licity is higher than 
is commonly assumed (e.g. lAronsll200ct ). Although internal 
dissipation in proto-magnetar winds is unlikely at t > t cn d, 
forced reconnection at large radii is a potential source of 
late-time emission in this case as well (see below). 

When the prompt emission ends a significant fraction 



10 Note that this problem does not arise in the standard internal 
shock model b ecause I" 1 , is assumed to v ary randomly throughout 
the burst (e.g. [Beloborodov et al. 2000), rather than to system- 
atically increase as predicted by the magnetar model. 



of the magnetar's initial rotational energy remains to be 
released in other forms. Since the beginning of the Swift 
mission, evidence has accumulated that GRB central en- 
gines are indeed active at late times, from minutes to J> 
hours following the burst. The X-ray afterglow in particu- 
lar shows a complex evolution, including a 'plateau' phase 
in the light curve which i s not predicted by the stan- 
dard forward shock model l|Nouse kl l200fj ; IWillingale et al.l 
120071 ). Superimposed on the smoother afte r glow are large 



ampl i tude X-ray flares (|Piro et all 120051; Burrows et al.l 



|2005| . I2OO7I ; IChincarini et all 120071 ; IChincarini et al.ll2010l ). 

which share many propert ies with the prompt GRB emis- 
sion |Margutti et al.ll2010l ) and als o appear to resul t from 
late-time central en gine activity |Lazzati fc Pernal 120071 ; 
iMargutti et~ai]|201fj| ). 

Although the magnetic dissipation or internal shocks re- 
sponsible for the prompt emission become ineffective when 
do is very large, spin-down luminosity can in principle power 
late-time emission in other ways. Indeed, a spin-down origin 
for the X-ray plateau is suggested by the 'plateau-like' evo- 
lution of the l ate-time wind powe r E(t) illustrated in Figures 
[2] and [5] (e.g. IZhang et al. 2006). Spin-down can in princi- 
ple power X-ray e mission either indirectly by refreshing the 
forwa rd shock (e.g. iGranot fc Kumarll2006l ; iDaH'Osso et al.l 
120101 ) or directly ('inter nally') by e.g. fo r ced reconnection 
at the forward shell (e.g lLvubarsWll2003l . 120051 ; iThompsonl 
120061 ; IZhang fc Yanll201ll) o r by upscattering forward shock 
photons IjPanaitescul 120081 ). Internal emission appears fa- 
vored in at least some cases due to the very steep decay 
observe d in the X-ray flux following the pla teau (e.g. GRB 
070110; iTroia et al.ll2007l ; iLvons et alj|20ld . hereafter L10; 
iRowlinson et alj|20ich l 

Figure [TH] shows a scatter plot of the wind power eval- 
uated at the beginning of the plateau-like high-cro phase 
^plateau = -E|t end as a function of the spin-down timescale 
r s | tcnd , calculated for models spanning the usual range of 
magnetar parameters (1 ms < Po < 5 ms; 3 x 10 14 G 
<; Bdip <; 3 x 10 16 G; x = 0,vr/2). Force-free spin down 
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Figure 18. Scatter plot of the wind power at the beginning of the plateau-like, high-cro phase -E p i a t C au = ^|t nd as a function of the 
spin-down timescale Ta|i end - Each point represents a model calculated within the range of initial spin periods 1 ms ^ Pq 5 ms and 
surface dipole fields 3 X 10 14 ^ B dip 3 X 10 16 G; results are shown for both mag netic obliquit i es x = aI1< i X = t/2. We also show 
for comparison the luminosities and end times t en( j x of the sample of plateaus from lLvons et all ((2010), which show a steep decline in 
flux at times t > t en< j x- Triangles and diamonds show the luminosities calculated assuming that the ratio between (observed) isotropic 
X-ray luminosity and wind power is equal to, or is a factor of 10 larger than, respectively, the gamma-ray beaming fraction (eq. [6]). 



is characterized by £ p i ateau oc Bdip^o.p and r = x B dil P o, P , 
such that r s oc -Bpiateau^Vp • The vertical scatter in Figure 
[l8]therefore results entirely from the distribution in 'initial' 
rotational periods Pq, p = P\t end following the GRB. 



L10 measure the isotropic X-ray luminosities Lx.iso 
and end times t cn a.x of the plateau phase for a subset of 
GRBs that show a steep decline in their X-ray flux at times 
£ > i C nd,x- In Figure [l8l we overplot i cn d,x and the luminosity 
from L10 Lx = ix.iso^x 1 corrected by a factor rjx = /b.xe^x 
that accounts for both the X-ray beaming fraction /b,x and 
the efficiency that spin-down power is converted into X-ray 
luminosity e ri x- We show two cases, in which rjx equals, or 
is a factor ~ 10 times larger than, the gamma-ray beaming 
fraction /b (which we estimate using equation ([6]) and the 
measured isotropic GRB energies). Note that because i en d,x 
is a lower limit on r s , figure [TS] shows that all of the plateaus 
measured by L10 are consistent with being powered by mag- 
netar spin-down for rjx <J 10/b- If t cn d is instead interpreted 



as the spin-down time itselfi"! our results indicate that ei- 
ther (1) the jet opening angle during the plateau phase is a 
few times larger than during the GRB itself, i.e. /b,x ^> fb 
and/or (2) the fraction of the spin-down power escaping 
through the jet and radiated in X-rays is 1. Although 
it is natural to expect that the radiative efficiency may be 
low when oo is very large at late times, too low of an effi- 
ciency may be inconsistent with afterglow energetics. It is 
also possible that a fraction of the late-time spin-down en- 
ergy is instead transferred to the supernova shock, although 
numerical simulations of the interaction of the wind with 
the st ar suggest this need no t be the case during the GRB 
itself l|Bucciantini et al.|[2009l ). 

Late-ti me magnetar ac t ivity could also produce X- 
ray flaring. iMargutti et al.l (|20 lOf ) find that the average 
flare luminosity decreases as Lflarc oc t~ a where a — 

11 As would be the case if spin-down triggers an abrupt end to 
the emission due to e.g. the delayed formation o f a black hole from 
a rot ationally-supported hyper-massive NS (e.g.rBaumgar te et al.l 
l2000l) . 
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2.7 (cf. lLazzati et al.ll2008l , iMargutti et al.ll2010l l. Although 
standard force- free spin-down predicts a = 2 at times 2> t s , 
steeper decays are inferred from the measured braking in- 
dices n of some pulsars (e.g. a = 4/(w — 1) ~ 2.42 fo r 
PSR J1846-0258 with n = 2.65; [Livingstone et all 120071 ). 
If prompt emission is indeed suppressed at late times by the 
high magnetization of the jet, periodic enhancements in the 
jet's mass-loading could temporarily 'revive' prompt-like in- 
ternal emission, resulting in flaring. Temporarily enhanced 
mass loss could result, for instance, from currents driven by 
a sudden rearrangement o f the magnetosphere, analogous t o 
Galactic magnetar flares (|Thompson &: Beloborodovll2005h . 
Indeed, X-ray flares could also be powered by the release of 
magnetic energy itself, which is ;> 10 49 — 10 50 ergs for typ- 
ical values of th e interior field strength B ~ 10 16 — 10 17 G. 
iGianniod rfiolph recently proposed searching for such 'super- 
flares' in nearby Galaxies, which could in principle be ob- 
served even long after the GRB, possibly in coincidence with 
a relic radio afterglow. 



6 DISCUSSION - A DIVERSITY OF 
PHENOMENA 

Magnetars may form with a variety of properties (and under 
a variety of conditions) which, in turn, manifests as a diver- 
sity of high energy phenomena. Figure [T5] summarizes the 
possible observable signatures of magnetar birth as a func- 
tion of the dipole field strength Bdi P and birth period Po. 
Although the plot shown is for an aligned rotator (x = 0) 
qualitatively similar results apply to the oblique case as well. 



6.1 Classical GRBs 

Magnetars in the upper left hand quadrant of Figure 1191 
produce 'classical GRBs' because (1) above the dotted lines 
the high energy emission is almost exclusively non-thermal 
because the relativistic jet dissipates its energy— through re- 
connection or shocks— above the photosphere beginning just 
after stellar break-out; (2) magnetars to the left of the dot- 
dashed line produce GRBs with energies E 1 > 10 50 ergs (see 
Fig. I10|) ; (3) magnetars in this regime produce outflow with 
average magnetization a avg ~ 10 2 — 10 3 , consistent with the 
inferred Lorentz factors of long GRBs (Figs. Illll2p . Note 
that the initial rotational energies of magnetars in this pa- 
rameter regime are > 3 x 10 51 ergs (Po < 3 ms), imply- 
ing that the requirements for a classical GRB and a hyper- 
energetic SN are remarkably similar (Fig. [T]). 



6.2 Very Luminous GRBs 

Magnetars in the extreme upper left corner of Figure [19] pro- 
duce classical GRBs with energies E 1 ~ 10 52 ergs which are 
comparable to the total rotational energy available (eq. [TJ). 
Evidence has recently grown for a class o f 'Very Luminous 
GRBs' (VLGRBs; e.g. lCenko et al.ll2010bl lah, which includes 
several Fermi bursts such as GRB 080 9 16C with an iso tropic 
energy £ 7 , iso w 8 x 10 54 ergs (e.g. lAbdo et al.ll2009l ). The 
observation that energetic Fermi bursts appear to be distin- 



guished by larger inferred Lorentz factorj 1 " 2 ! Y > 10 3 than is 
estimated for more typical GRBs is consistent with the cor- 
relation predicted by the magnetar model between the aver- 
age GRB luminosity and jet magnetization <T aV g (maximum 
Lorentz factor), as shown in Figure 1171 Many extremely 
energ etic GRBs, such as GRB 990123 (e .g. Kulkarn i et al.l 
1 19991) and 080319B (|Bloom et all 120091 1. are also distin- 
guished by bright optical emission coincident with the GRB. 
The synchrotron emission predicted by the magnetic dissi- 
pation model at optical-UV wavelengths contributes an es- 
pecially large fraction of the total radiated energy in bursts 
with large magnetization er avg (see the late-time spectra in 
Fig. [15]). 

At present, the properties of VLGRBs appear consistent 
with resulting from magnetars with extreme, but physically 
reasonable, properties. However, measurements of the total 
energy in relativistic ejecta E to t — E 1 + E^ (where Ek is 
the kinetic energy) could constrain— or even rule out— the 
magnetar model as the central engine if E to t were found 
to exceed the maximum rotational energy ~ -Erot(Po ~ 
1 ms) ~ 3 x 10 52 ergs. Although efforts are presently under 
way to determine E to t for a sample of well-studied bursts 
jCenko et al.ll2010bl lah. the results of these studies are hin- 
dered at present by simplifying assumptions in the after- 
glow modeling and jet structure , which may lead to system- 
atic overestimates in £ /tot (e.g. Izhang fc MacFadvenl [20091 ; 
Ivan Eerten et al"1l2010T l. Nevertheless, VLGRBs provide im- 
portant probes of the most extreme central engine proper- 
ties. 



6.3 Low Luminosity GRBs 

Magnetars to the right of the dot-dashed line in Figure 
1191 produce GRBs with energies <J 10 50 ergs which may 
contribute to t he class of so-called ' low luminosity GRBs ' 



(LLGRBs; e.g. iBloom et al.l |2003bl: ISoderberg et al.l 120041 : 

" 2Q7C 



ICobb et al.ll2006l ; lLiang et al.ll2007l ; lKaneko et alj2007t ). L~L 



GRBs are distinguished from classical GRBs by their lower 
energies, simple gamma-ray light curves (generally a single 
pulse), longer durations, and higher local ra tes fe.g. lCowardl 
120051 ; iLe fc Permed 120071 ; iLiang et al.ll2007h . Because large 
angular momentum is probably rare in core collapse su- 
pernovae, LLGRB-producing magnetars with weaker fields 
and/or slower rotation may indeed be formed more com- 
monly than the magnetars responsible for classical GRBs. 



6.4 Thermal-Rich GRBs and X-Ray Flashes 

Magnetars below the dotted lines in Figure Q1J] produce jets 
that dissipate a significant fraction of their energy under 
optically thick conditions after breaking through the star 
(i.e. they pass through Stage IVa described in JJJ) and pro- 
duce jets with lower Lorentz factors than classical GRBs, 
i.e. cr avg < 10 2 . We speculate that proto-magnetars in this 
regime may produce X-ray-rich GRBs or X-ray Flashes 
fXRFs; iHeise et alJlioOll ; iMazzali et al.ll2006l ') because they 
are accompanied by lower-frequency, quasi-thermal emission 



12 Note, however, that the lower limit constraints on T derived 
for Fermi bursts become weaker if the ~ GeV and y MeV photons 
originate from different radii (e.g. IZou et aljfeoiof) . 
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Figure 19. Regimes of high energy phenomena produced by magnetar birth in core collapse supernovae, as a function of the magnetic 
dipole field strength Bjip and initial rotation period Pg, calculated for an aligned rotator (\ = 0). 



with an energy comparable to, or somewhat lower than, the 
non-thermal GRB emission itself (Fig. [9]). Although XRFs 
share many properties with long GRBs, such as an associ- 
ation with mass i ve st a r form ation (e.g. iBloom et alJ 120031 : 
ISoderberg etafl |2004 120071 ). they may be distinguished 
from GRBs by their ability to cou ple a significant energ y 
to highly relativistic material fe.g. ISoderberg et alJ I2004T 1. 
This is consistent with the fact that magnetars in the lower 
portions of Figure [19] indeed radiate a smaller fraction of 
their total energy during the GRB (as compared to the 
radiatively-inefncient high-cro phase; Q than magnetars in 
the classical GRB regime. 



6.5 Choked Jets and Very Luminous Supernovae 

Magnetars in the lower right hand corner of Figure [19] 
produce jets with peak isotropic luminosities <J 10 48 ergs 
s -1 . Low power jet may be unstable (|Bucciantini et aT] 
2009) or take longer to propa gate through the star than 
the duration of the GRB (e.g. iMeszaros fc Waxmanl l200ll ; 
IWaxman fe Meszarosl [200I ) . Magnetars in this regime may 
thus produce 'choked' jets with little direct electromagnetic 
radiation (although they could still be a source of high en- 



ergy cosmic rays or neutrinos; e.g. [Waxmanl [l995l ; IVietril 
1 19951 ; lAndo fc Beacomll2005l ; iMurase et al.ll2009r i. 



A number of core-collapse SNe have been recently dis- 
cove red that are unusually bright and/or optically-energetic 
Ofek et all [2007]; ISmith et all 1200?]; ISmith et af 



Quimbv et aU 120071: iRest et al.ll2009l ; iGal-Yam et al 
Quimbv et al.l 20091 ) . Proposed explanations for these 



2008; 



2007; 



events, collectively known as very luminous SNe (VLSNe), 
include pair-instability SNe (|Barkat et al.l Il967l ); interac- 
tion of the supernova shock wit h dense circumstell a r ma- 
terial (e.g. IGal-Yam et ail 1200?] ; ISmith et all |2007t 120081 ; 
lMetzgerll2010l ); and the injection of late-t ime rotational en- 
ergy from a rapidly- spinning magnetar l|Kasen fc Bildstenl 
120101 ; IWooslevI 120101) . In order to energize the supernova 
ejecta on t he ~ days-weeks timesca les relevant for power- 
ing VLSNe, iKasen fc Bildstenl |2010T ) conclude that a mag- 
netar with Bdip ~ 5 x 10 14 G must possess an initial rota- 
tion period Po ~ 2 — 20 ms. This nominally places VLSNe- 
producing magnetars in the 'choked jet' regime. We note, 
however, that in order to explain VLSNe, the initially Poynt- 
ing flux-dominated magnetar wind must thermalize its en- 
ergy behind the SN sh ock, instead of escaping in a jet 
jBucciantini et al"1l2009h . which might still be able to prop- 
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agate through the star on the longer timescales of relevance 
for VLSNe. 



6.6 Galactic Magnetars 

If known Galactic magnetars were born with magnetic fields 
similar to their current observed strengths Bdip ~ 10 14 — 10 15 
G (e.g. iKouveliotou et ail 1 19981 ') and as fast rotators, then 
Figure [19] suggests that their formation was accompanied 
by a thermal-rich GRB/XRF or choked jet, depending on 
their initial rotational period. Slower rotation, correspond- 
ing to a choked jet, may be likely in the majority of cases 
because Gal actic magnetars are formed in ~ 10% of core 
collapse SN (| Woods fc Thompson! lioOrj ) . yet only a small 
fraction of en velope-stripped SN are accompanied by rela- 
tivistic ejecta (|Soderberg et al]|2006h . Furthermore, the SN 
remnants of Galactic magnetars do not show evidence for 
hyper-energe tic SN explosions (e.g . IVink fc Kuiper] 120061 ; 
see, however. iHorvath fc AllerfeoiOT l. 



6.7 Accretion-Induced Collapse 

This paper has focused on the core collapse of massive stars, 
but magnetars may also form via the a ccretion-induced col- 
lapse (AIC ) of white dwarfs (WD, e.g. iNomoto et~ai1ll979l ; 
lUsovl fl992T l. Although AIC is probably intrinsic ally rarer 
than standard core collapse (e.g. iFrver et al.lll999h . millisec- 
ond magnetars may be a more common byproduct of AIC 
because the WD is spun up considerably as it accretes up 
to the Chandrasekhar mass. 

A distinguishing characteristic of AIC is the lack of a 
massive overlying stellar envelope. However, AIC does not 
produce a vacuum around the magnetar. A small quantity 
of mass ~ 10~ 3 — 10~ 1 M [ n is eject ed during the supernova 
explosion itself l|Dessart et al.f2006r i and in the early, mildly- 
relativistic phase of the neutrino wind (Stage II). If the col- 
lapsing white dwarf furthermore has sufficient angular mo- 
mentum, an accretion disk form s around the neutron star 
|Michej|l987l ; lDessart et alfeOOfj ). As this disk accretes onto 
the NS on a timescale <C 1 s, outflows from the disk powered 
by nuclear recombination eject J> 10~ 2 Mq in Nickel-rich 
material (|Metzger et al.ll2009l )Fl 

Because the proto-magnetar is surrounded by a modest 
'sheath' of material, its relativistic wind from the magnetar 
may be collimated into a bipolar jet, analogous to the stan- 
dard core collapse case. Because of the lower inertia of this 
surrounding mass, however, collimation may be less effective 
and the opening angle of any 'jet'-like structure may be con- 
siderably larger. If this speculation is correct, it would imply 
a larger beaming fraction /b, lower isotropic luminosity, and 
softer spectral peak (e.g. eq. [11]) than in the core collapse 
case. Perhaps equally important, the fact that the jet is no 
longer required to escape the star in order to produce high 
energy emission may 'select' for magnetars with lower fields 
and/or slower rotation (and, hence, lower spin-down lumi- 
nosities, lower Lorentz factors, and softer spectra) than in 
the core collapse case. 

13 This implies that although AIC is not accompanied by a 
bright supernova, it may produce a dimmer ~ day-long tra nsient 
jMetzger et al.ll2009l; IFrver et al1l2009l; iDarbha etal]|2010h . 



One of the biggest mysteries associated with short- 
duration GRBs is that J> 1/4 are followed by a 'tail' of emis- 
sion (usually soft X-rays) starting ~ 10 seconds after the 
30 — 100 seconds ([Norris fc Bonnell 



2006 




2009 


: 



Gehrels et al.ll200d ; ICal-Yam et al]|2006l : iPerlev et al 



Norris et alJl201Cft . Although the large inferred ener- 



gies and durations o f the tails are difficul t to explain in NS 
merger models (e.g. iMetzger et alj 120101 ). their properties 
are similar to the prompt em i ssion expected from magnetar 
birth via AIC. lMetzger et alj l|2008h proposed an AIC model 
for 'short GRBs with extended emission', in which the short 
GRB is powered by the accretion of the disk onto the NS as 
described above and the subsequent 'tail' is powered by the 
(wider-angle) proto-magnetar wind. This model is consistent 
with the host galaxy demographics, and the la ck of a bright 



supernova, associated wi t h short GRBs (e.p 
20061 ; iBerger et all 120051; iHiorth et al.1 12005 



Bloom et al 



Perlev et al 



20091 : lBergerll2010] ; lFong et alj|2010h . We note that an anal- 
ogous model invoking a long-lived magnetar remnant that 
survives a NS-NS merger could also in principle explain the 
late-time X-ray activity. Such a possibility is supported by 
the re cent discovery of a « 2Mq pulsar bv lDemorest et all 
(2010), which demonstrates that the high density equation 
of state is relatively stiff. 



7 CONCLUSIONS 

In this paper we take the first steps towards developing the 
millisecond proto-magnetar model into a quantitative theory 
for gamma-ray bursts. Using detailed evolutionary models of 
magnetar spin-down, we explore a wide range of magnetar 
properties and calculate the prompt emission predicted by 
magnetic dissipation and internal shock models. Although 
the picture we construct may not be accurate in all details, 
it serves as a 'proof of principle' that the basic concepts 
can be constructed into a self-consistent model. Our work 
also provides a baseline for future improvements, as will be 
necessitated in particular by advances in our understanding 
of the origin of prompt GRB emission. 

Several theoretical uncertainties remain that should be 
addressed with future work. These include a more detailed 
understanding of the effects of rotation and convection on 
the cooling evolution of the proto-neutron star, and the ef- 
fects of strong magnetic fields on the neutrino-driven mass 
loss rate. Although most of our results are at least qual- 
itatively robust to these uncertainties, predictions for the 
GRB duration (and how it correlates with other observ- 
ables) is in particular sensitive to the time of neutrino 
transparency. The mass loss rate from the proto-NS (and, 
hence, the wind magnetization) also depends on fraction of 
the magnetosphere open to outflows, which depends on the 
poorly-understood sources of dissipation near the Y-point. 
Future studies would also be aided by a more detailed un- 
derstanding of the dependence of the jet properties (e.g. 
break-out time and opening angle) on the properties of the 
proto-magnetar and the stellar progenitor. The source of 
the rapid rotation and strong magnetic fields required to 
produce millisecond magnetars also remains a major un- 
certainty. However, we note that black hole models place 
similar, if not more extreme, constraints on the progenitor 
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rotation and the large-scal e magnetic field of the central en- 
gine fe.g. lMcKinnevll2006l ). 

Our primary conclusion is that a surprisingly large frac- 
tion of GRB properties can be explained by the magnetar 
model. These include: 

• Energy. Magnetars with properties in the 'classical 
GRB' regime in Figure \W\ radiate E~, ~ 10 50 - 10 52 
ergs during the GRB phase, consistent with the beaming- 
corrected gamma-ray energies inferred from afterglow mod- 
eling. Magnetars with stronger (weaker) magnetic fields 
and/or shorter (longer) initial periods may produce very lu- 
minous(low luminosity) GRBs. 

• Lorentz Factor. Magnetars in the 'classical GRB' 
regime produce jets with average and instantaneous mag- 
netizations do > 10 2 — 10 3 (Fig. 1 1 1 p which are remarkably 
similar to the typical Lorentz factors inferred from GRB 
observations (cf. Fig. I12p . The baryon loading of the jet is 
not fine-tuned or put in by hand, but instead results natu- 
rally from the physics of neutrino heating above the proto- 
magnetar surface. This is contrast to black hole models, for 
which current predictions for Y depend on the uncertain 
rate at which baryons diffuse into an otherwi se clean jet 
(e.g. iLevinson fc Eichlerll2003l ; iMcKinnevI [20051 ) . The mag- 
netar model predicts that ao (and probably V) increases 
monotonically with time during the burst. Among other 
things, this implies that any thermal emission present will be 
strongest at early times and will decrease in relative strength 
as the outflow becomes cleaner with time (Fig. I15[) . 

• Duration. The GRB begins once the jet breaks out of 
the star and becomes optically thin at the internal shock 
or magnetic dissipation radius. The GRB ends once the jet 
magnetization increases sufficiently that jet acceleration and 
dissipation become ineffective. Because the latter generally 
occurs when the NS becomes transparent to neutrinos at 
t = t„-thin ~ 10 - 100 s (Fig. the magnetar model 
naturally explains the typical durations of long GRBs. 

• Steep Decay Phase. The abrupt onset of the high-<To 
transition at t « i„-thin (Fig. [2} explains why GRB prompt 
emission decreases rapidly a fter the prompt emissi on ends 
(the 'steep decay' phase; e.g. iTagliaferri et al.ll2005l ). 

• Association with Hypernova. It is natural to as- 
sociate energetic, MHD-powered supernovae with magne- 
tar birth. If the magnetar model is correct, all long GRBs 
formed from the core collapse of massive stars should be 
accompanied by an energetic (and possibly hyper-energetic) 
supernova. Magnetars formed via AIC, by contrast, may pro- 
duce long GRBs not accompanied by a bright SN. This is a 
promising explanation for the ~ 100 second X-ray tails ob- 
served following some short GRBs ( i|6.7[) and explains why 
they resemble long GRBs in many of their properties. 

• High Lorentz FactorsoEnergetic Bursts. The 
magnetar model predicts a positive correlation (with signif- 
icant scatter) between the (energy-weighted) average mag- 
netization (Javg of the jet and the (beaming-corrected) GRB 
luminosity /energy (Fig. I17[) . This is consistent with the 
fact that energetic Fermi bursts appear to have the largest 
Lorentz factors. 

• High Radiative Efficiency. Both magnetic dissipa- 
tion and internal shocks may occur in proto-magnetar winds, 
resulting in the prompt high-energy emission. Both mod- 
els predict maximum radiative efficiencies e r ~ 30 — 50%, 



consistent with the high values of e r inferred from after- 
glow modeling (e.g . Panaitescu fc Kumaijl200l1 ; IZhang et al.1 
l2007l ; lFan fc Pirarj|2006t ). 

• Amati-Yonetoku Relation. Our spectral modelling 
favors magnetic dissipation over internal shocks as the 
prompt emission mechanism, in part because magnetic 
dissipation predicts a relatively constant spectral energy 
peak i?p C ak as a function of time (Fig. I16[) . Strong inter- 
nal shocks may be suppressed by the residual magnetiza- 
tion of the ejecta or if the toroidal fi eld geometry is not 
conducive to particle acceleration (e.g. ISironi fc Spitkovskvl 
[2010)). In combination with the predicted o" a v g — L 1 corre- 
lation fFig. I17[) . the magnetic dissipation model reproduces 
both the slope and normalization of the observed Amati- 
Yonetoku correlations. 

• Late-Time Emission. Although we expect that 
prompt internal emission becomes ineffective when ao be- 
comes very large at late times, the plateau X-ray afterglow 
phase may also be powered by magnetar spin-down, as pro- 
posed by previous authors and suggested by Figure [2] The 
predicted correlation between the plateau luminosity and 
duration (Fig. I18[) is c onsistent with the sample of 'inter- 
nal' plateaus studied bv lLvons et al.1 (|2010l ) . Late-time X-ray 
flaring may be powered by residual rotational or magnetic 
energy. 
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APPENDIX A: PROTO-MAGNETAR WIND 
PROPERTIES 

In this appendix we describe how to calculate the power E 
an d mass loss rate M f t]A2[) of proto-magnetar winds. 



Al Energy Loss Rate E 

The winds from millisecond proto-magnetars are acceler- 
ated primarily by magnetic forces rather than by thermal 
pressure. At large radii the wind power can be divided into 
components of kinetic energy and magnetic Poynting flux, 

viz. E — -Ekin + S ma g. 

The kinetic luminosity of the wind is -Ekin = {Too — 
l)Mc 2 , where Too = (1 — f^/c 2 ) -1 / 2 and Voa are the asymp- 
totic Lorentz factor and velocity of the outflow, respectively. 
Nonrelativistic outflows have magnetization ao < 1 and 

1/3 

reach an asymptotic speed Voo ~ c<r , resulting in a kinetic 
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Table 1. Properties of Proto-Magnetar Winds 
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( a ) Surface dipole magnetic field strength following NS contraction. C^Spin period if NS were to contract to final radius with fixed 
angular momentum [Actual initial spin period at t=0]. ( c ) Magnetic obliquity (see Fig. [3jl . W'Stretch' correction applied to neutrino 
luminosities and energies to account for the effects of rotation (see eq. IA110 , ( e 'Wind magnetization when the jet breaks out of the 
stellar surface at t = i bo = 10 s (eq. \5\ ; Fig. [6j . Energy- weighted average magnetization between jet break-out and the end of the 
prompt emission, t DO < t < *cnd (Fig- [Til l. Energy- weighted average Lorentz factor of the bulk shell produced by internal shocks 
(Fig. |12|). (^Duration of the prompt GRB emission Tgo = t cn d — *thin is (Fig. [l3j) . ( l 'Time after core bounce when the prompt GRB 
emission ends, defined as the point when the 'saturation radius' r mag (eq. 1101 ) exceeds the internal shock radius r s h (eq. IB3h . This 

transition generally occurs simultaneous with the transition of the proto-NS to neutrino transparency (see Fig. [JJ. ( J ) Maximum 
'thermal' energy produced by the jet (eq. [8]). ^Maximum GRB energy, denned as the rotational energy released in the time interval 
min[t^ , tthin.is] < * < *end ( see Fig. 1100 . (^Wind power at t ~ t cn d. ( m ) Dipole spin-down tim escale at t = t on( j. ^ Calculated using the 

neutrino cooling calculations of Hiidcpo hl et al.l ||2010|) . 



luminosity -Bkin oc (1/2) Mt& l|Lamers fc Cassinellil Il999h . 

1/3 

Relativistic outflows have ao > 1 and achieve Too « a 
near the fast magnetoson ic surface, beyond which accel- 
eration effectively ceases (|Goldreich fc Julian! fl97oh . This 
weak 1/3 power embodies the classical problem that (uncon- 
fined, time-sta tionary) high-op winds ac celerate inefficiently 
in ideal MHD ijKennel fc Coroniti|[l983 : see gj]for further 
discussion). 

The magnetosphere is completely open outside of the 
Alfven radius Ra ■ The Poynting flux -E mag is thus related to 
the open magnetic flux cj>, and hence to the magnetization 
ctq (eq. [2] ) , via the relationship 



T-\E x B 



2<j> 2 fl 



x rdQ 



2 / 2 2^ 

3 (v^B^r ) 



Ra 



3Ra 



Mc 2 <r 2/3 , 
Ale 2 CTQ, 



(70 > 1 



,(A1) 



where E = — (v /c)xB is the electric field and the factor 2/3 
accounts for the angular integration. The equalities in the 
second line follow because (1) the outflow co-rotates with 



the star, such that ~ fir out to radii ~ Ra, where v<j, 
is the toroidal velocity; and (2) near Ra the poloidal field 
begins to bend back appreciably due to the fluid's inertia. 
The toroidal magnetic field strength £?0 thus becomes com- 
parable to the poloidal field ~ £? r at r ~ Ra, such that 
4> = B T r 2 \n A ~ B < f > \n A R 2 A (the equality is exact in the case 
of force- free winds) . In writing the third and fourth lines we 
have made use of the fact that 

Vao/Q. ~ Rl^a^ 3 , (70 <C 1 

Rl = c/n (70 > 1 



Ra = 



(A2) 



The open magnetic flux of a rotating dipole with a 
surface magnetic field strength Bdi p is given by 



j 'open B dip Rn 



(A3) 



where 



(1/2tt) / dep sm6 d9 

Jo Jo 



1 — COS 



2 J 



(1 + sin x) 



2R Y 
(A4) 
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is the fraction of the NS surface threaded by open field. Here 
# opcn ~ 2sin -1 [(J? ns /JiY) 1/ ' 2 ] is the opening angle of the po- 
lar cap, corrected by a factor (1 + sin 2 x) 1 ^ 2 t o account for 
the larger cap size of an obli que rotator (e.g. ICheng et al.l 
l200d ; lBai fc Spitkovskvll201(]| ); R ns is the NS radius; R Y is 
the 'Y' point radius where the close zone ends in the mag- 
netic equatorial plane; and the second equality holds in the 
small-cap limit # opon <C 1 (Ry 2> -Rns). See Figure [3] for an 
illustration of the relevant geometry. 

Just after core bounce, thermal pressure may dominate 
above the NS surface and the entire magnetosphere may 
open into a 'split-monopole' configuration with / opC n ~ 1. As 
the NS contracts, cools, and spins up, however, its magnetic 
field is amplified and magnetic pressure eventually comes 
to dominate. This produces a 'closed' or 'dead' zone at low 
magnetic latitudes from which a steady-state wind cannot 
escape (i.e. Ry > Rns)- In the limit of a force- free wind 
(cro ^ 1) the radius of the Y-point likely extends close to the 
radius of the light cylinder, but in general Ry is ^ Rl for 
less m agnetized (finite-ao) winds. Following iMetzger et al.l 
( 2007 1, we assume that Ry/Rl = min[0.3oo 15 , 1] for R Y > 
Rns, based on an empirica l fit to the axisymmetr ic rela- 
tivistic MHD simulations of iBucciantini et all |2006h . which 
span the non-relativistic to relativistic transition. The val- 
ues of Ry that we adopt are si milar to those w e estim ate 
by applying the toy model of iMestel fc Spruitl (| 19871 ') to 
the proto-magnetar context. Determining the detailed time- 
dependence of -Ry will, however, ultimately require incorpo- 
rating a self-consistent, physical model for the resistivity in 
the magnetosphere and equatorial current sheet. 

Combining our results, the total wind power is given by 

E — Bmag + -E"kin 

_ f McVf, a < 1 
\ (2/3)Mc 2 a , (T > 1 ' 

(A5) 

where the magnetization (eq. [2]) can now be written 

d2 r>4 q2 r2 
^dip ^ns 1 6 J open 

°~o — ^— 

Mc 3 

B 2 ip R ns » 4 (l + sin 2 X ) (Ry\~ 2 (A6) 
fl Y »-R„ 5 4Mc 5 V Rl / 

Note that in the nonrelativistic case the kinetic and mag- 
netic contributions to the total power are similar (-E mag = 
2-Ekin), while in the relativistic case the outflow is Poynting 
dominated since -Bmag/^kin ~ o" 2 / 3 ^ !• 

A2 Mass Loss Rate M 

Mass loss from the proto-NS results from neutrino heating 
in the atmosphere just above the NS surface. The dominant 
heating and cooling processes are the charged-current reac- 
tions 

u e + n o e~ + p and v e + p -o- e + + n (A7) 

For unmagnetized winds, the mass lo ss rate is well- 
appro ximated by the analytic expression ijQian fc Wooslevl 
1 19961 ) 

/ j V 5/3 

M„ = 5 x 1O" 5 M s" 1 ——-Z r x 

w \ 10 52 ergss -1 / 



( e v V'V3 ( M ns V 2 ( R ns \ s / 3 5/3 

{ lOMeV J { OTq ) { iOknJ (1 + £cs) " (A8) 

Although both electron neutrinos and antineutrinos con- 
tribute to the heating, for simplicity we combine their con- 
tributions into a single product of the neutrino luminosity 
L v and mean energy e„, defined by 

L„e 2 = Z/„Je„J 2 + Z^Je^J 2 , (A9) 

where the |...| represent an appropriate average over the neu- 
trino absorption cross sections. The normalization adopted 
in equation (|A8[l inclu des both this averaging and a general 
relativistic correction l|Thompson et al.ll200lh . The parame- 
ter 

•-^(^(isErfe)" 1 (Mo) 

is a correction <C 1 for t he additional heating d ue to inelastic 
electron scattering fsee lQian fc Wooslevlll996l . eq. 50). 

In mos t calcu lations we use L v (t),t v (t), and R ns (t) from 
iPons et al l \ 19991 ). hereafter P99, who calculate the delep- 
toni zation and cooling evolutio n of non-rotating proto-NSs 
(cf. iBurrows fc Lattimerl [l986l ). Examples of L u (t),e u (t), 
and R ns (t) are shown in Figure IXT1 for different NS masses. 
Note that for t J> 1 s, L v and e v decrease relatively gradually 
as a power-law until a time t v -tbXn ~ 10 — 60 seconds, after 
which L v and e„ plummet as the proto-NS becomes trans- 
parent to neutrinos. As we show in |J3] tv-thm determines 
the GRB duration in the proto-magnetar model. 

Since M v depends sensitively on L v and e„ we briefly 
discuss the limitations and the uncertainties in the calcu- 
lations of P99. First, although portions of the proto-NS 
are c onvectively unstable du r ing its early co oling evolu- 
tion (IBurrows fc Frvxelll [l993l ; iKeil et"allll996h . convective 
transport is not accounted for by P99. The primary effect 
of convection is to increase the cooling rate and hence to 
speed up the temporal evolution of the neutrino luminosity 
(L. Roberts, private communication). P99 find that the rate 
at which L v and e„ decrease at late times, and hence the 
precise value of t„_thin, also depends sensitively on the high 
density equation of state, which is uncertain. In order to 
explore the sensitivity of our results to uncertainties in L v 
and e v , we also calculate models using neutrino luminosities 
a nd energies from the re cent proto-NS cooling calculations 
of lHudepohl et all l|2010h . hereafter H10 (L. Roberts, private 
communica tion), which follow a successful electron-capture 
supernova l|Kitaura et al.l [20061 ). This calculation, which in- 
cludes improvements in the neutrino opacities over previous 
work, is shown for comparison in Figure IA1I The primary 
difference between the cooling curves of P99 and H10 is the 
significantly faster late-time evolution found by H10. 

Finally, neither P99 or H10 include the effects of mag- 
netic fields or rotation, yet this paper focuses on proto- 
magnetars rotating at a significant fraction of their break- 
up speed. Rapid rotation decreases the interior temper- 
ature of the NS, which slows its cooling evolution. Us- 
ing one-dimensi o nal r otating core collapse calculations, 
iThompson et all l|2005h find that L v and e„ are reduced 
by factors of ~ 0.5 and ~ 0.8, respectively in their fastest 
rotating model at t « 600 ms following core bounce, com- 
pared to an otherwise equivalent nonrotating case. Ideally 
the effects of rotation on L v and e v should be calculated self- 
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Figure Al. —L eft Panel: Weight ed neutrino luminosity L v (defined in eq. |A9| ) as a function of time after core bounce. The first three 
models are from lPons et al. calculated for NS masses M ns = 1.2Mq (solid line), 1.4Mq (dotted line), and 2.0Mq (dashed line). 

The last model is from iHudepohl et al] l|2010h (M ns = 1.27Mq; dot-dashed line). Right Panel: Mean electron anti-neutrino energy e v 
(left axis) and NS radius (right axis) as a function of time after core bounce, for the same models shown in the left panel (the radius 
from H10 is not shown). 



consistently. Lacking such a model, however, we account for 
rotational effects qualitatively by introducing a 'stretch' pa- 
rameter 77 s , which modifies the cooling evolution from the 
non-rotating case (fi = 0) as follows: 



L v — > L v \n=oWa 



e w |n=o>7< 



-1/4 



(All) 



where a value r/ a ~ few is motivated by the calculations of 
iThompson et all ([2005) for millisecond rotators. Note that 
this simple parametrization preserves the total energy radi- 
ated in electron neutrinos and increases the time of neutrino 
transparency tp-tun oc r/ s . Although we expect rj a to be an 
increasing function of fl, in our calculations we fix r) s = 3 
for lack of a predictive model. We also neglect differences in 
the neutrin o radiation field with latitude caused by rapid ro- 
tation (e.g. iBrandt et al]|2010» . which if properly included 
would impart the total wind mass loss rate with an addi- 
tional dependence on the magnetic obliquity \ 

A strong magnetic field modifies M from the standard 
expression in equation (|A8[) in three ways. First, M is re- 
duced by a factor / opC n (eq. |A4| ~) since only the open frac- 
tion of the surface contributes to the outflow. Second, M is 
enhanced by a factor /cent due to centrifugal 'slinging.' This 
occurs when rotation is sufficiently rapid and the magnetic 
field is sufficiently strong that centr ifugal forces increase th e 
scale height in the heating region ^Thompson et all 12004 



By fitting the numerical results of iMetzger et al 



we find that the maximum centrifugal enhancement to M 
(obtained in the strong-field limit of strict co-rotation de- 
scribed below) is well-approximated by the functional form 
(for P > 1 ms) 



/ccnt.max ^ exp[(P c /P)^] 

for a value /3 — 1.5, where 



(A12) 



P c ~ 2.1 sin a 



llOknJ {TaM^J mS ' (A13) 

where a « max[6 , opcn /2, \\ is a typical angle from the rota- 
tional axis sampled by the open zone. The normalization we 



ad opt for P c is d e termin ed by fitting the numerical results 
of IMetzger et all (|2008T ). which were calculated for equa- 
torial field lines (a = tt/2). The scaling of P c with mass, 
radius, and a, however, are chose n based on the theoret i- 
cal expectation that P c oc Rj_/c s (|Thompson et all |2004| ) . 
where R± ~ R ns sin a and c s are the centrifugal 'lever arm' 
and sound speed in the gain region, respectively. The latter 
is proportional to the NS escape speed oc (M/_R ns ) 1//2 (see 
iQian fc Wooslevl[l99rj . eq. 45). 

Although /ccnt.max is the maximum enhancement of M, 
it obtains only if the magnetic field is sufficiently strong 
that the outflow co-rotates with the star to a location out- 
side the sonic radius R s . This requires Ra > R B = (t^t) 1 '' 3 
llLamers fc Cassinellil fl999l ) . Using the numerical results of 
IMetzger et alj (|2008l ). we find that a satisfactory interpola- 
tion of the mass loss enhancement between the centrifugal 
(-Ra S> -Rs) and non-centrifugal (Ra Rs) regimes is given 
by /cent = /ccnt,max(l-exp[-R A /-Rs]) + exp[-_RA/Rs]. Note 
that relativistic outflows are necessarily in the centrifugal 
regime because for oo ^> 1, Ra ~ -Rl = c/fi > -R s (eq. |A2j ). 

Finally, a strong magnetic field changes M v (eq. |A8| 1 by 
altering the neutrino heating and cooling r a tes in the proto- 
NS a tmosphere (e.g. iDuan fc Qianl |2004| ; iRiquel me et al.l 
2005). The most important effect is that the electrons 
and positrons participating in the charged-particle reac- 
tions (eq. IA7I ) are restricted into discrete Landau levels 



i|Duan fc Qianll2004h . In this paper we neglect these effects 
because we estimate that the corrections to M„ are rela- 
tively minor for surface field strengths -Bdip ^3x 10 16 G0 
However, more detailed future work should address the de- 
pendence of M v on Bdip. 



14 Note that M v depends most sensitively on the heating and 
cooling rates in the gain region, which is typically located a 
few kilometers above the NS surface, where the magnetic field 
strength B oc r~ 3 is a factor of a few weaker than the surface 
field strength. 
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Once the NS becomes transparent to neutrinos at late 
times, L„ and e„ decrease rapidly (Fig. IA1|I and the mass 
loss rate decreases abruptly. Neutrino heating only deter- 
mines M so long as the magnetosphere is sufficiently dense 
that vacuum electric fields do not develop. This assumption 
breaks down, however, on ce M decreases to near the critical 
iGoldreich fc Julian! l| 19691 ) (GJ69) mass loss rate. However, 
because in actual pulsar winds M exceeds the GJ69 value, 
we instead assume that the minimum mass loss rate in the 
pair-dominated regime is given by a multiple of GJ69 rate, 



/'- 



+m e (7/e) 

Bdip 



3 x 10" 
P 



s 



(A14) 



(^-+\ ( ^di P \ (p_y ( R ™ \ 

{ 10 6 Al0 15 G/ V msj VlOknJ 

where I = 47t_Rl 7 ?gj|r l c, m c and e are the electron mass 
and charge, and r/GJ ~ (flB/2-Kc) is the GJ69 charge den- 
sity, evaluated at the light cylinder. The multiplicity fi |_ of 

positron/electrons produced by magnetospheric acceleration 
is un certain, especially in the case of magnetars l| Thompson! 
2008). Lacking a predictive model, in our calculations we 
fix the multiplicity at a value = 10 6 which is consis- 

tent with estimates based on detai led synchrotron emissio n 
models of pulsar wind nebulae (e.g. lBucciantini et al.|[201Cf ). 
Although the late-time wind magnetization depends sensi- 
tively on the multiplicity (o"o|t»t„_ thin oc l//i_ + ), most of 
our conclusions regarding late-time emission (ij5]) are insen- 
sitive to this choice. 

To summarize, the mass loss rate is given by 



M 



where 



en j 

■iV/y j open j 'cent ; 

Mgj, 



(X + fopen/2) < 7T/2 
(X + ^open/2) £ Tv/2 

M ^ M G j 



(A15) 



/cent = /ccnt,max(l - exp[-i? A /-Rs] ) + CXp [- R A / R s ] (A16) 

and 



Ra „ ( ape 3 \ 1/3 
R s \ GMQ. j 



(A17) 



APPENDIX B: INTERNAL SHOCK EMISSION 

In this Appendix we derive the luminosity and syn- 
chrotron spectral peak of internal shock emission from proto- 
magnetar winds as described in H4.4I 

In a small interval dt centered around the time ij , the jet 



releases a mass dM s = M- s \t- 3 dt with energy dE s 
Lorentz factor Fj L . The collision between this fast shell and 
the bulk shell occurs at the time f s h and radius 7? s h defined 
by the condition 

rt s h 

R sh = j p s dt = (U-*j)/9j, 



Ej\ tj dt and 



(Bl) 



-2x1/2 



ifTj.r, 



where # = (1 - P. 
times) and Tj ^> T s (valid at late times) one finds that 

r 2 

L s 
f2 

I 

2f 2 



tsh — tb. 

rj»r a 
Rsh ^ 2 



2r 2 tj 



S> 1 (usually valid at all 



(B2) 



^ ( l - 2ff ) 



(B3) 



Because both r s and tj increase with time, Rsh also in- 
creases as the GRB proceeds. Neglecting cosmological di- 
lation, emission from the shock is received by a distant ob- 
server at the time 



to 



Rsh/c + tbo — tj+tbo, 

r J »r B 



(B4) 



where in the second equality we have used equations (|B2|) 
and (|B3[) . This shows explicitly how relativistic effects con- 
spire to produce emission on timescales matching those set 
by the central engine. When the fast shell becomes fully in- 
corporated into the bulk, the bulk Lorentz factor increases 
by an amount dr s determined by conservation of momen- 
tum, 



2M S 



(B5) 



resulting in the release of thermal energy 



dE sh = M s F s c 2 + dMjFjc 2 
dM ]C 2 ( r 2 



(dMj+M s )(r s + dr s 



2r, 



(B6) 



If a fraction e c of dE^ is imparted to the electrons, which 
then radiate it away on the expansion timescale, the total 
radiative efficiency is given by 



dE s i x e & 
dMTjC 2 



r 2 r 

1 + — - 2— 

+ r 2 r 

J 



(B7) 



At early times Pj > T B and e r <l, but at later times when 



P 3> r s , e r reaches a maximum value 



3 /2. 



We now consider how the jet and shell interact in 
greater detail in order to evaluate the peak energy of the re- 
sulting synchrotron emission. In the frame of the bulk shell 
(hereafter denoted by a tilde), the jet velocity and Lorentz 
factor are given by 



ft 



r s 2 /r 2 



i + r 2 /r 2 



f j = (i - 



-1/2 



r;»r B 



ij/2r B 



(B8) 



Note that unlike in the standard internal shock scenario, the 
relative Lorentz factor between the shells Fj is 3> 1 at late 
times when Fj 3> F s . This allows for high radiative efficiency. 

If the shocked gas is relativistically hot, the (rest frame) 
post-shoc k number and energy dens ities are given, respec- 
tively, by (jBlandford fc McKedll976l ) 



n sh = (4fj + 3)nj; e sh = (fj - l)n sh m p c 2 
where 



(B9) 



(BIO) 



is the pre-shock density and Lj,i so = Ljf^ 1 . 

If a fraction e e and £(, of the post-shock energy is par- 
titioned into electron kinetic and magnetic energy, respec- 
tively, the resulting (electron) random Lorentz factor and 
magnetic field strength are given by 



e e m p (rj - 1) 
7e = : h 1 



B 



e e m p Fj 

.ISO 



rj»r s 2( e m e T s 

1/2 



(Bll) 



1/2 r 1/2 
,21/2rfc3/2 t 



(B12) 
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where C, e is the fraction of electrons accelerated. 

The peak synchrotron photon energy as seen by the 
observer is then 

m e c rj>r B 
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